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Abstract 

The  Doppler-broadened  two-photon  absorption  (TPA)  cross-section  for  the 
52S'i/2  — *  52D5/2  transition  in  Rb  is  measured  using  direct  absorption  methods.  The 
selection  rule  \AF\  <  2  applied  to  both  isotopes  yields  17  transitions  in  3  Doppler 
limited  lines.  A  detailed  model  of  the  intensity  profile  was  also  developed  to  account 
for  a  focused  Gaussian  beam  (with  an  M2  value  of  1.09)  propagating  through  a  two- 
photon  absorption  medium.  A  peak  absorbance  of  24%  was  observed  for  an  intensity 
of  6.28  ^2  at  the  focus,  a  Rb  density  of  4.6  x  1015  cm~ 3,  and  a  path  length  of 
15  cm.  Alkali  concentrations  from  1.61  —  8.52  x  1015  cm~3  were  monitored  in  the 
far  wing  of  the  D2  line.  Extracting  the  hyperfine-broadened  TPA  cross-section  from 
87  test  configurations,  while  varying  the  pump  power,  alkali  concentration  and  focal 
length,  yielded  an  error-weighted  average  of  6.75  x  10~21^^  with  a  standard  deviation 
of  3.61  x  10” 21  •  This  cross-section  is  sufficient  for  a  pulsed  dye  laser  to  bleach 

the  pump  transition  in  the  Two-Photon  Pumped  Alkali  Laser  (TPAL)  that  lases  at 
420  nm  and  5.2  jim. 

Optically  pumped  atomic  rubidium  lasers  pumped  in  the  blue  satellite  of  the  D2 
line  from  the  ground  Rb-Ar  or  Rb-Kr  collision  pair  to  the  dissociative  R2E^2  state 
produce  laser  emission  at  780.2  nm.  Lasing  is  achieved  for  pump  wavelengths  of  752.3 
to  greater  than  760  nm  for  the  Rb-Ar  system  and  757.1  —  760.4  nm  for  the  Rb-Kr 
system.  Slope  efficiencies  increase  with  both  Rb  and  Ar  concentrations  and  exceed 
0.25%  using  a  heat  pipe  configuration.  The  gain  is  very  high  with  photon  build-up 
times  of  1  —  3.7  ns.  Laser  induced  heating  and  subsequent  condensation  of  alkali  vapor 
in  the  heat  pipe  configuration  currently  limits  operation  to  less  than  2500  Torr. 
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OPTICALLY  PUMPED  ATOMIC  RUBIDIUM  LASERS:  TWO-PHOTON  AND 
EXCIPLEX  EXCITATION  MECHANISMS 

I.  Introduction 

Ever  since  the  effective  demonstration  of  the  Diode-Pnmped  Alkali  Laser  (DPAL) 
by  Krnpke  et  al.  in  2003  [65],  the  laser  community  has  focused  much  of  their  attention 
on  investigating  this  field  in  the  hopes  of  developing  a  non-chemical  high  power  laser. 
The  alkali  laser  is  very  promising  for  many  reasons.  First  of  all,  the  DPAL  is  not  a 
chemical  laser,  so  it  does  not  require  the  storage  or  the  replenishment  of  highly  toxic 
chemicals.  Instead,  the  DPAL  lasing  medium  is  the  alkali  vapor  contained  neatly 
and  compactly  inside  a  glass  cell  the  size  of  a  spool  of  thread.  This  alkali  vapor  is 
also  reusable  as  the  same  atom  is  cycled  through  the  process  100  million  times  every 
second.  The  DPAL  system  is  also  not  a  solid  state  laser,  so  it  does  not  have  the 
thermal  issues  attributed  to  such  a  device  as  pump  powers  are  increased.  In  fact, 
there  is  minimal  heat  generated  in  the  entire  process  since  the  only  heat  generated 
originates  from  the  relaxing  of  the  atom  from  the  n2P3/2  state  (where  n=3,4,5  &  6 
corresponding  to  Na,  K,  Rb  &  Cs  respectively)  down  to  the  n2P!/2  state.  This  gap 
represents  238  cm-1  to  554  cm~ 1  for  rubidium  (Rb)  and  cesium  (Cs)  respectively  and 
these  are  both  on  the  order  of  magnitude  of  kT  at  room  temperature. 

The  mechanics  of  the  DPAL  system  are  also  quite  simple.  Using  the  cesium  atom 
as  an  example,  the  pump  laser  operating  at  852  nm  excites  the  ground  state  (62S'i/2) 
atom  to  the  62P3/2  state.  Collisional  mixing  then  relaxes  the  atom  to  the  62Pi/2  state. 
Population  inversion  then  occurs  between  the  62Pi/2  state  and  the  ground  state  where 
lasing  occurs  at  895  nm.  This  particular  wavelength  is  in  the  near-infrared  band.  For 
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a  more  detailed  discussion  of  this  process  please  refer  to  the  papers  prepared  by 
Krupke  et  al.  [65,  66]. 

More  recently,  higher  order  variants  to  the  DPAL  have  been  reported  to  lase  back 
to  the  ground  state.  These  variants  may  draw  population  away  from  the  regular 
DPAL  but  also  increase  the  number  of  wavelengths  at  which  DPAL-like  lasers  have 
access  to.  This  research  reports  findings  on  two  variants,  the  Exciplex  Pumped  Alkali 
Laser  (XPAL)  and  the  Two-Photon  Pumped  Alkali  Laser  (TPAL). 

The  TPAL  is  pumped  through  a  virtual  state,  halfway  between  the  alkali  ground 
state  and  the  target  pump  level.  As  the  alkali  atom  encounters  a  photon  at  the 
frequency  of  the  virtual  state,  the  atom  simultaneously  absorbs  two  identical  photons 
moving  population  to  the  real  pump  state.  In  2008,  Sulham  et  al.  demonstrated  that 
a  Cs  TPAL  system  will  simultaneously  yield  a  ’’blue”  and  a  rnid-IR  lasing  wavelength 
[16,  40,  55]  in  a  cascade  effect  down  to  the  ground  state.  The  mechanism  for  achieving 
population  inversion  in  this  upper  state  is  more  involved  than  the  originally  designed 
DPAL  system,  but  the  applications  for  such  a  device  are  desirable1.  This  research 
reports  the  absolute  two-photon  absorption  (TPA)  cross-section  for  the  52Si/2  — > 
52D5/2  transition  in  a  Rb  vapor,  in  an  effort  to  understand  pump  intensity  thresholds 
for  a  possible  Rb  TPAL.  Although  the  TPA  cross-section  for  this  transition  in  Rb  was 
first  measured  in  1983  by  Zapka  et  al.,  there  has  been  subsequent  groups  reporting 
values  that  differ  by  as  much  as  a  factor  of  100.  To  restore  confidence  in  the  true 
value,  this  research  varied  alkali  concentration,  the  beam  waist  and  pump  power  to 
generate  87  different  test  configurations.  The  resultant  average  number,  weighted 
against  the  error  of  each  test  case,  is  reported  and  compared. 

The  XPAL  system  is  pumped  in  the  blue  satellite  wing  of  the  D2  line  and  lases 
either  on  the  D2  or  Dx  line  [34,  35,  36,  37,  38].  A  pulsed  XPAL  was  first  demonstrated 

^uch  applications  are  under  water  communications,  low  diffraction  HEL,  and  IR  counter  mea¬ 
sures  at  multiple  wavelengths 
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in  2009  by  pumping  Cs-Ar  excimers  in  a  mixture  of  500  Torr  of  argon  and  100  Torr 
of  ethane.  The  excited  excimcr  then  populates  the  upper  laser  level  (2Pi/2)  in  atomic 
cesium  by  decay  [35].  Lasing  occurs  on  the  Di  line  at  894.3  nm  after  spin-orbit 
relaxation,  aided  by  the  ethane.  The  excimer  step  creates  a  broad  band  of  pump 
wavelegths  (~  5  nm)  and  removes  the  requirement  for  hydrocarbons  [34,  35,  36, 
37,  38].  The  four-level  Cs-Ar  system  exhibits  a  threshold  at  130  —  250  J  in  a  4  ns 
pulse  with  an  efficiency  of  about  1%  and  has  been  scaled  in  pump  rate  to  2-3  times 
threshold  [4,  36,  37].  This  research  extends  the  four-level  XPAL  work  by  lasing  in 
Rb-Ar  and  Rb-Kr  mixtures  and  reports  slope  efficiencies  and  laser  output  powers. 
This  research  was  the  first  to  use  a  heat  pipe  configuration  for  the  XPAL,  and  was 
the  first  to  explore  rare-gas  pressures  of  up  to  2500  Torr. 
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II.  Previous  Work 


The  motivation  for  the  Rb  TPAL  originated  from  an  experiment  performed  in  the 
summer  of  2009  by  a  team  at  the  Air  Force  Institute  of  Technology  (AFIT).  During 
the  course  of  tuning  a  Spectra  Physics  Nd:YAG  pumped  Sirah  dye  laser1  to  various 
excited  states  in  atomic  rubidium  in  an  energy  pooling  experiment,  a  bright  blue  dot 
was  noticed  parallel  to  the  pump  beam.  It  was  later  realized  that  this  occurred  when 
the  laser  was  tuned  slightly  off  resonance  near  the  52Ps/2  state  in  Rb.  Due  to  the 
nature  of  the  experiment,  the  rubidium  cell  was  only  being  interrogated  by  a  pump 
laser,  so  no  resonator  existed  around  the  cell.  It  was  reported  that  ASE  occurred  at 
about  420  nm  and  in  the  Mid-IR  in  a  range  from  3.85  -  5.2  [im. 

The  motivation  for  the  Rb  XPAL  came  from  the  joint  work  performed  by  the 
DPAL  group  at  AFIT  and  the  XPAL  group  at  the  University  of  Illinois,  Urbana- 
Champaigne  (UIUC)  [33].  The  same  Spectra  Physics  Nd:YAG  pumped  Sirah  dye 
laser  above  was  used  to  pump  a  17  cm  long  pyrex  cell  with  a  Rb-Ar  mixture  inside. 
Slope  efficiencies  were  investigated  maintaining  constant  temperature  and  pressure. 

2.1  The  One-Photon  Pump  Alkali  Lasers 

Krupke’s  Diode-Pumped  Alkali  Laser  (DPAL). 

In  2003,  William  Krupke  patented  the  design  of  the  first  diode-pumped  alkali  laser 
system  [65].  Although  this  was  a  first- of- a-kind  laser  system,  the  alkali  vapors  were 
well  known  by  spectroscopists.  As  early  as  1958,  Schawlow  and  Townes  [61]  used  the 
potassium  vapor  to  push  the  maser  technology  into  the  optical  and  IR  wavelengths. 
In  1962,  two  years  after  the  establishment  of  the  laser  [51],  Rabinowitz  et  al.  was 
the  first  to  demonstrate  an  alkali  vapor  as  an  actual  gain  medium  in  his  continuous 
1This  laser  delivered  up  to  100  mJ  in  10  ns  at  10  Hz. 
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optically  pumped  cesium  laser  [54],  This  concept  did  not  gain  a  strong  footing  in  the 
laser  community  since  little  can  be  found  in  the  literature  afterwards.  This  is  probably 
due  to  the  lack  of  efficient  pump  diodes  at  the  time.  For  the  next  four  decades,  the 
alkali  vapors  moved  away  from  being  the  gain  medium  to  becoming  the  the  work 
horses  in  areas  such  as  metrology  and  spectroscopy.  High-power  pulsed  and  cw  lasers 
radiated  sealed  glass  cells  containing  alkali  vapors  while  the  side  fluorescence  was 
collected  and  analyzed  to  determine  important  atomic  and  molecular  properties.  The 
cesium  atom  in  particular  is  probably  the  most  characterized  atom  next  to  hydrogen 
as  the  unit  of  time,  the  second,  is  defined  in  terms  of  9,192,631,770  periods  of  the 
radiation  corresponding  to  the  transition  between  the  two  hyperhne  levels  of  the 
ground  state  of  cesium  [13].  As  a  result,  the  second  has  become  the  unit  with  the 
smallest  uncertainty.  The  affinity  towards  the  alkali  in  the  scientific  community  is 
due  primarily  to  their  atomic  structure.  They  possess  one  valence  electron  making 
it  hydrogen-like,  such  as  having  a  smaller  than  usual  ionization  energy.  Compared 
to  other  atoms  it  doesn’t  take  a  significant  amount  of  energy  to  excite  the  atom. 
Therefore,  the  chemistry,  spectroscopy  and  the  kinetics  of  the  alkalis  were  all  a  fairly 
mature  concept  by  the  time  William  Krupke  developed  the  DPAL. 


Figure  1.  This  is  the  three-level  diagram  for  the  DPAL  system  using  the  rubidium 
vapor  as  the  gain  medium  [66]. 
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Krupke’s  DPAL  functions  in  the  following  manner.  All  alkali  vapors  have  a  ground 
atomic  state  of  2Si/2,  with  the  next  two  excited  states  being  2P\/2  and  2P3/2.  There  are 
two  electric-dipole  allowed  transitions  in  this  three-level  system.  These  transitions, 
the  2S'i/2-2Pi/2  and  the  2*Si/2-2P3/2,  are  more  commonly  called  the  Di  and  the  D2  lines 
respectively  (see  Figure  1).  The  D2  transition  is  the  pump  transition.  So,  the  pump 
laser  for  a  rubidium  cell  would  be  tuned  to  the  780  nm  wavelength  in  order  to  move 
population  from  the  ground  state  to  the  52P3/2  level.  The  natural  lifetimes  for  the 
52P3/2  and  the  52Pi/2  states  are  26.2  ns  and  27.7  ns  respectively  [64],  The  population 
would  normally  bottle-neck  in  the  upper  state  because  the  spin-orbit  relaxation  rate 
is  slow  compared  to  the  pump  rate.  To  compensate  for  this,  a  buffer  gas  is  added  on 
the  order  of  1  atmosphere  to  increase  the  relaxation  rate  from  52P^/2  down  to  52Pi/2. 
Finally,  with  population  inversion  achieved  between  the  2P\/2  state  and  the  ground 
state,  lasing  in  the  DPAL  occurs.  At  this  point  the  process  is  able  to  start  all  over 
again  for  the  alkali  atom,  since  the  pump  laser  is  still  radiating  the  cell.  For  a  more 
in-depth  study  of  the  DPAL  laser  please  refer  to  the  following  articles  by  Krupke, 
Beach,  Page  and  Zhdanov  [11,  12,  56,  58,  65,  66]. 

One  of  the  acknowledged  challenges  to  the  DPAL  system  is  the  narrowness  of  the 
spectral  width  of  the  D2  transition  (approx.  10  GHz  at  1  atmosphere).  This  becomes 
a  problem  when  developing  a  high  power  laser,  because  the  spectral  width  from  bars 
and  stacks  in  a  semiconductor  pump  laser  at  high  powers  is  about  100  times  larger 
than  that.  One  way  engineers  have  tried  to  mitigate  this  is  to  broaden  the  spectral 
width  of  D2  by  adding  inert  gases  to  the  alkali  vapor.  This  broadening  comes  at  no 
small  cost  since  enough  inert  gas  must  be  added  to  increase  the  pressure  up  to  25-50 
atmospheres. 
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2.2  Two-Photon  Absorption  Spectroscopy 


The  theory  of  two-photon  absorption  using  a  virtual  intermediate  state  is  not 
new  to  the  scientific  community.  In  fact,  the  theory  is  just  about  as  old  as  the 
theory  of  quantum  mechanics  itself.  In  1931,  Maria  Goppert-Mayer  published  an 
article  in  the  Annalen  Der  Physik  pioneering  the  framework  of  the  theory  of  two- 
photon  absorption  [25]  from  the  then  new  field  of  quantum  mechanics.  This  theory 
remained  unproven  for  three  decades  waiting  for  the  invention  of  the  laser2  in  1960 
by  Maimen[51].  The  theory  of  two-photon  absorption  was  finally  confirmed  in  an 
experiment  using  an  europium-doped  calcium  fluoride  crystal  by  Kaiser  et  al.  in 
1961  [44],  and  then  in  a  cesium  vapor  by  Abella  in  1962  [8].  For  the  next  several 
years  many  articles  were  published  exploring  this  new  field  of  spectroscopy.  Much 
of  the  research  focused  on  multiphoton  ionization,  whether  it  was  of  hydrogen  and 
rare-gas  atoms  [14]  or  of  the  alkali  atoms  [15].  There  was  also  some  work  done  in 
1965  and  1967  investigating  the  infrared  emission  during  the  two-photon  excitation  of 
the  62S'i/2  state  in  potassium  vapors  [60,  68].  Recall  that  the  group  at  AFIT  in  2009 
concluded  that  there  were  IR  emission  in  their  two-photon  excitation  experiment  with 
a  rubidium  vapor.  Another  more  attractive  attribute  of  the  two-photon  absorption 
spectroscopy  is  the  ability  to  directly  access  states  that  are  forbidden  in  one  photon 
transitions  [22], 

The  1970’s  saw  the  emergence  of  two-photon  Doppler-free  spectroscopy.  This 
allowed  some  states  to  be  directly  accessible  which  were  once  inaccessible  due  to  parity 
incompatibility  [29] .  This  also  allowed  the  hyperfine  structure  of  high-lying  states  to 
be  more  accurately  resolved  and  consequently  the  hyperfine  splitting  to  be  determined 
with  greater  fidelity.  One  such  study  was  done  by  a  group  at  Harvard  University  in 

2 Although  Gordon  Gould  published  the  term  LASER  in  1959[26],  the  laser  was  also  referred  to 
as  an  optical  maser  at  this  time. 
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1974  [47].  They  studied  the  32S'i/2-52S'i/2  transition  in  a  sodium  vapor  and  determined 
the  hyperfine  splitting  for  the  52<Si/2  state  to  be  101  ±  15  M Hz.  With  the  invention  of 
the  Ti:Al203  (Ti-sapphire)  laser  in  1982  [52]  the  field  of  spectroscopy  reached  higher 
levels  of  fidelity.  In  contrast  to  the  dye  lasers,  the  Ti-sapphire  laser  had  a  tuning  range 
that  covered  three  regions  of  three  different  dyes,  and  was  able  to  remain  at  high  power 
for  longer  periods  of  time  without  degradation  in  beam  quality.  In  addition  to  this, 
the  Ti-sapphire  laser  reached  higher  output  power  and  a  narrower  linewidth.  Needless 
to  say,  this  opened  the  door  to  better  spectroscopy.  In  1990,  Adams  et  al.  used  a  Ti- 
sapphire  ring  laser  to  perform  high  resolution  spectroscopy  of  the  52S'i/2  <—  52-D5/23/2 
transition  [9]  in  rubidium.  Adams  et  al.  knew  that  the  Doppler  broadened  linewidth 


85Rb87Rb  87Rb  85Rb 

0  5  10  GHz 

Figure  2.  The  rubidium  52Si/2  — >  52D3/2  transition  obtained  by  Doppler-free  two-photon 
absorption  spectroscopy  with  a  Tksaphire  laser  [9]. 

of  the  52S'i/2  —>  5 2D3/2  transition  at  room  temperature  would  be  about  0.5  GHz. 
Therefore  they  chose  to  use  Doppler-free  two-photon  absorption  spectroscopy.  Their 
results  are  shown  in  Figure  2.  Also  in  1990,  Garreau  et  al.  [32]  used  Doppler-free 
two-photon  spectroscopy  on  atomic  hydrogen  and  deuterium  to  gain  a  more  accurate 
value  of  the  Rydberg  constant,  Roc3.  They  were  able  to  achieve  observed  frequencies 
with  an  accuracy  of  1  part  in  109.  In  2000,  the  Rydberg  constant  was  again  put 


3The  Rydberg  constant  is  the  most  accurately  known  physical  constant  known  to  man. 


to  the  test  when  Niering  et  ah  measured  the  absolute  frequency  of  the  hydrogen 
12S'i/2  <—  22S'i/2  two-photon  transition  with  an  accuracy  of  1.8  parts  in  1014  [50]. 

With  the  more  recent  development  of  laser  cooling  and  trapping  techniques,  two- 
photon  spectroscopy  has  been  able  to  join  magneto-optical  trap  systems  to  measure 
more  precise  two-photon  ionization  cross-sections  [23].  In  2009  a  non-linear  optics 
team  from  the  University  of  California  expanded  the  generalized  Kramers- Heisenberg 
expressions  to  include  two-photon  absorption  [53].  Also  in  the  late  2000’s,  two-photon 
absorption  has  been  seen  as  an  emerging  player  in  the  realm  of  quantum  computing 
[27].  All  of  these  examples  indicate  a  need  for  recognizing  the  presence  of  and  the 
mechanism  behind  two-photon  absorption. 

2.3  The  Two-Photon  Pump  Alkali  Lasers  (TPAL) 

Sulham  et  ah  reported  ASE  from  one  or  both  5 2P  states  when  pumping  the 
virtual  states  that  correspond  to  the  52D3/2, 52H5/2  and  the  72Si/2  states  (Figure  3). 
This  group  was  also  the  Erst  to  construct  a  resonator  based  on  the  two-photon  pump. 
However,  three  unexpected  results  occurred  after  constructing  such  a  cavity.  For  one, 
the  threshold  decreased  as  the  concentration  of  the  alkali  increased.  Slope  efficiencies 
tend  to  increase  with  alkali  concentration,  but  it  is  anomalous  for  the  threshold  to 
decrease.  Secondly,  it  was  shown  that  the  resonator  actually  degraded  the  power  of 
the  ASE  line.  This  too  is  counterintuitive  since  the  gain  of  the  laser  should  increase 
as  the  optical  path  length  increases.  Thirdly,  it  was  observed  that  the  power  out 
decreased  as  the  concentration  of  He  increased.  This  is  unusual  because  the  spectral 
width  of  the  two-photon  pumped  n2Si/2  —  n2D5/2,3/2  transition  is  6  times  smaller  than 
the  spectral  width  of  the  pump  laser,  which  in  this  case,  is  16  GHz  delivered  by  a 
Spectra  Physics  Nd:YAG  pumped  Sirah  dye  laser  [16].  The  He  (or  other  rare  gas) 
broadens  the  spectral  width  of  the  transition  to  better  match  the  spectral  width  of 
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the  pump  laser. 


26,311.44 


Figure  3.  The  Grotrian  diagram  for  rubidium.  The  solid  arrows  are  the  two-photon 
pump  transitions.  The  dotted  arrows  are  the  amplified  stimulated  IR  emission  lines. 
The  dashed  arrows  are  the  amplified  stimulated  ’’blue”  optical  emission  lines  [16]. 


2.4  The  Two-Photon  Absorption  Cross-Section 

The  two-photon  absorption  (TPA)  cross-section,  a2<t>,  is  an  important  quantity  to 
know  in  order  to  calculate  the  pump  intensity  threshold  for  a  potential  two-photon 
absorption  laser  system.  For  the  transition  in  this  research  there  are  two  simultaneous 
laser  transitions  back  to  the  ground  state,  see  Figure  4,  that  can  be  taken  advantage 
of  by  a  Two-Photon  Pumped  Alkali  Laser  (TPAL). 

Referring  to  Figure  4,  the  following  rate  equations  can  be  derived  to  attain  an 
approximation  of  the  threshold  condition. 


dn3 

dt 

dn2 

dt 


&2  <j> 


ni  -  r IRn3 


T  irti3  —  A21n2 


(1) 

(2) 


where  is  the  branching  ratio  for  the  first  laser  transition  and  A2 1  is  the  sponta- 
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Figure  4.  The  three  laser  levels  for  the  Two-Photon  Pumped  Alkali  Laser  (TPAL). 

neous  emission  for  the  upper  ’’blue”  laser  level.  Applying  a  steady  state  condition  to 
Equation  1,  we  can  define  an  approximation  to  the  laser  threshold  condition  to  be, 


Assuming  that  we  have  no  losses  in  the  system,  a  — >  0,  we  can  define  the  popu¬ 
lation  ratio  in  Equation  3  to  be  equal  to  the  degeneracy  ratio,  or  2.  Equation  3  now 
shows  the  importance  of  knowing  the  TPA  cross-section  as  a  first  step  in  developing 
a  TPAL. 

There  have  been  a  few  reports  on  the  value  of  the  two-photon  absorption  cross- 
section  of  the  transition  52S'i/2  — »  5 2D5/2  for  rubidium.  Table  1  displays  the  values 
of  the  TPA  cross-section  for  the  525i/2  — >  52_D5/2  transition  from  the  literature. 
It  is  clear  from  this  table  that  the  value  of  the  two-photon  absorption  cross-section 
is  far  from  being  well  established.  It  is  one  goal  of  this  research  to  establish  a  firm 
understanding  of  the  true  value  by  performing  a  comprehensive  study  of  this  transition 
in  Rb. 

The  first  value  reported  was  in  1983  by  Zapka  et  al.  His  group  used  Doppler- 
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Table  1.  Two-photon  absorption  cross-sections  measured 
for  the  52/Si/2  — ►  52D5/2  Transition  in  Rb. 


(  cm 4  \ 

Reference 

Year 

Type 

Isotope 

4xl0_2U 

Zapka  et  al.  [67] 

1983 

Experimental 

87Rb 

1.2xl0-18 

Collins  et  al.  [19] 

1993 

Experimental 

85Rb 

0.57xl0-18 

Marinescu  et  al.  [49] 

1994 

Theoretical 

85Rb 

6.5x  10-20 

Saha  et  al.  [43] 

2011 

Theoretical 

85Rb 

5  xlO"20 

Saha  et  al.  [43] 

2011 

Experimental 

85Rb 

free,  frequency-modulation  (FM)  spectroscopy  to  detect  the  weak  absorption  in  the 
52S'i/2  — >  52D,5/2  transition  in  an  isotopically  enriched  Rb  sample  (number  densities 
85 Rb  :8'  Rb  «  1  :  6)  [67].  A  droplet  of  this  isotropically  enriched  Rb  was  inserted  into 
a  glass  cell,  which  was  contained  in  an  oven  heated  by  coaxial  wires  and  monitored 
by  a  thermocouple.  Referring  to  Figure  5,  a  vertically-polarized  single-mode  cw 
pump  laser  (of  about  120  mW )  was  set  at  frequency  u0  fa  12851  cm-1  and  traversed 
a  polarizing  beam  splitter  and  a  15  cm  focal  length  lens  before  entering  the  test 
cell.  After  recollimating,  the  beam  passed  a  phase  modulator  that  is  driven  by  an 
rf  oscillator  at  a  frequency  vrf  between  100  and  600  MHz.  After  going  through  a 
quarter- wave  plate  the  beam  is  reflected  by  an  etalon,  which  acts  as  a  perfect  mirror  at 
the  two-photon  pump  frequency  and  becomes  the  probe  beam  back  through  the  test 
cell.  FM  spectroscopy  was  performed  by  scanning  the  laser  frequency  vq  and  keeping 
urf  constant.  The  absolute  absorption  of  this  transition  was  determined  by  calibrating 
the  FM-absorption-shaped  signal  amplitude  on  a  known  reference  absorption  signal. 
This  was  done  by  tuning  to  one  of  the  adjacent  resonant  frequencies  of  the  etalon. 
This  procedure  then  yielded  an  absolute  absorption  of  2.5%  for  the  52S'i/2(A  =  2)  — > 
52A>/2(F  =  4)  transition  in  Sl  Rb.  A  temperature  of  182  °C  amounted  to  a  number 
density  of  about  3.3  x  1014  cm"3.  The  carrier  beam  was  estimated  to  be  about 
60  mW,  with  a  beam  waist  of  70  gm.  Under  these  conditions  an  approximate  value 
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°f  cr 2(b  =  4  x  10~20  was  determined  and  reported.  Besides  the  Doppler-free 
techniques  of  counter-propagating  beams,  the  experimental  set-up  in  Zapka’s  work 
most  closely  resembled  the  set-up  in  our  present  work. 


Polarizing 
beam  splitter 


Figure  5.  The  apparatus  used  by  Zapka  et  al.  [67]. 

The  second  value  wasn’t  reported  until  10  years  later  by  Collins  et  al.  This 
group  used  Doppler  free,  differential  absorption  spectroscopy  (DAS)  with  an  external- 
grating-cavity  cw  diode  laser,  set  at  100  to  tune  through  the  52Si/2(F  =  2)  — > 
52D5/2(F  =  0,1,  2,  3, 4)  transitions  of  85 Rb  in  a  2.54  cm  long  Pyrex  cell  containing 
a  number  density  of  6.8  x  10~18  cm~ 3  [19].  The  absorption  signal  is  detected  by  an 
electronic  noise-canceling  circuit  that  cancels  up  to  60  dB  of  common-mode  noise. 
The  laser  frequency  was  swept  with  a  triangle  function  at  400  Hz.  The  resulting 
absorption  signal  has  a  FWHM  of  43  MHz ,  which  is  about  22  MHz  larger  than  the 
range  over  which  the  hyperhne  transition  is  distributed.  It  was  reported  that  the 
added  linewidth  was  contributed  by  a  slight  Doppler  broadening  (due  to  the  coun- 
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terpropagating  beams  coming  in  at  an  angle  of  0.95°)  and  laser  frequency  jitter.  Key 
to  determining  the  TPA  cross-section  in  this  group  was  calculating  the  effective  two- 
photon  line  strength,  S ^2\  In  reaching  this  number,  Collins  determines  the  integrated 
absorption  coefficient  and  the  total  atomic  excitation  rate.  The  integrated  absorp¬ 
tion  coefficient  is  found  by  integrating  the  absorption  curve  determined  by  fitting  five 
Gaussian  functions  (corresponding  to  the  five  hyperfine  lines)  with  the  same  FWHM 
to  the  absorption  feature.  Accuracy  of  the  final  number  in  this  method  becomes  de¬ 
pendent  on  the  knowledge  of  the  factors  involved  for  the  FWHM  of  the  absorption 
feature.  Maintaining  a  temperature  at  143  °C  (i.e.  N  =  6.8  x  1013  cm~ 3)  Collins  et 
al.  reported  a  number  that  was  two  orders  of  magnitude  larger  than  Zapka  et  al.  ten 
years  earlier,  =  1-2  x  10~18  ^ . 

One  year  later,  Marinescu  et  ah  reports  a  theoretical  value  based  on  modeling 
reduced  absorption  rates  for  this  TPA  transition  [49].  The  number  that  this  group 
reported  was  <j2(p  =  0.57  x  10~18  g^.  Unfortunately,  the  method  Marinescu  chose 
makes  it  extremely  difficult  to  compare  with  this  work’s  reported  value. 

Most  recently  in  2011,  a  TPA  cross-section  was  measured  by  Saha  et  ah  using  low 
powers  (1  mW )  and  low  densities  in  a  hollow-core  photonic-band-gap  fiber  (PBGFs) 
filled  with  Rb  vapor  [43].  A  9  —  cm  long  and  6  fim  diameter  hollow-core  PBGF 
was  mounted  in  a  vacuum  chamber  connected  to  an  Rb  source  containing  naturally 
occurring  Rb.  The  temperature  of  the  chamber  and  finger  were  maintained  at  85  °C 
and  55  °C  respectively.  After  bake-out  the  background  pressure  was  measured  at 
10  8  Torr.  An  external  cavity  diode  laser  beam  was  split  50/50  to  counterpropagate 
through  the  fiber  of  equal  power  and  identical  polarization.  To  generate  the  desired 
OD  100)  a  highly  off-resonant  20  mW  desorption  beam  at  805  nm  was  also 
coupled  to  the  fiber  with  an  orthogonal  polarization.  They  measured  the  TPA  as 
a  function  of  the  intensity  of  light  using  an  acouto-optic  modulator  (AOM)  to  vary 
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the  intensity  of  the  two  counter-propagating  beams  simultaneously  using  a  triangle 
wave  at  1  kHz.  The  emitted  blue  photons  are  collected  and  detected  from  the  top 
of  the  fiber.  For  a  direct  measure  of  the  TPA  from  the  778.1  nm  pump  beams,  they 
performed  a  lock-in  detection  by  using  another  AOM,  to  modulate  one  of  the  two 
counter-propagating  beams.  A  weak  reflection  of  the  forward  propagating  beam  is 
collected  by  a  photodiode.  A  1%  absorption  was  determined  by  detecting  30  nW 
of  output  blue  flourescence  at  420.3  nm  with  a  1  rnW  diode  pump  laser.  With 
this  method  a  direct  absorption  value  was  calculated  to  be  5  x  10~20  for  the 
52S'i/2  — >  52D5/2  in  85Rb.  This  value,  although  an  order  of  magnitude  different  from 
Marinescu,  came  to  within  20%  of  Zapka’s  reported  value. 

It  is  clear  from  the  literature  that  an  accepted  value  for  the  TPA  cross-section 
for  this  particular  transition  in  Rb  has  not  been  agreed  upon.  As  stated  earlier, 
developing  a  TPAL  on  this  transition,  requires  knowledge  of  the  TPA  cross-section. 

The  Excimer  Pumped  Alkali  Laser  (XPAL). 

Referring  to  Figure  6,  the  rubidium  and  argon  combine  to  form  the  Rb-Ar  excimer 
in  it’s  ground  state,  X2E%2.  The  pump  laser  excites  the  excimer  to  the  blue  wing  of 
the  repulsive  P2E(y2  state  which  then  quickly  dissociates  into  a  ground  state  argon 
atom  and  the  excited  52P3/2  state  of  atomic  rubidium.  The  remaining  process  mimics 
that  of  a  regular  DPAL,  where  the  ethane  aides  in  the  spin  orbit  relaxing  of  the  52P3/2 
state  to  the  upper  laser  lever,  the  52Pi/2  state.  The  fifth  level  is  the  atomic  ground 
state  of  Rb  completing  the  XPAL  process. 

In  the  latter  half  of  2008,  Readle  et  al.  from  the  University  of  Illinois,  devised 
a  method  to  circumvent  the  spectral  width  problem  mentioned  above  in  the  DPAL 
system.  Instead  of  moving  population  from  the  ground  state  of  an  atomic  alkali  atom, 
they  pumped  the  ground  state  of  an  excimer  in  the  blue  satellite  (837  nm)  to  the 
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Internuclear  Distance,  R  (A) 

Figure  6.  The  first  four  Rb-Ar  surfaces  calculated  by  Blank  et  al.  [46]. 

unstable  B 2X%2  state  that  decayed  into  the  atomic  2P3/2  state.  The  first  experiments 
using  this  technique  had  500  Torr  of  Argon  and  100  Torr  of  ethane  in  Cs-Ar-C2H6 
mixtures  [35,  36,  38].  It  was  reported  that  this  pumping  mechanism  produced  broad 
spectral  linewidths  5  nm ),  and  slope  efficiencies  approaching  10%  with  absorbed 
pump  energy  thresholds  of  ~  40  /iJ .  This  was  first  done  using  a  pulsed  dye  laser, 
driven  by  the  second  harmonic  of  an  Nd:YAG  system  generating  4  —  5  ns  FWHM 
pulses.  With  pump  energies  as  large  as  6  mJ ,  the  dye  laser  has  a  pulse  rep  rate  and 
nominal  linewidth  of  10  Hz  and  ~  0.2  cm~2  respectively.  The  alkali  cell  used  in  their 
experiment  was  made  of  quartz,  had  brewster  angles  and  had  dimensions  of  6.0  cm 
in  length  and  2.5  cm  in  diameter.  This  XPAL  system  used  five  levels  to  generate  the 
894.3  nm  (the  Cs  Dj  line)  lasing  line.  The  first  two  of  the  five  levels  are  the  ground 
and  excited  states  of  the  Cs-Ar  excimer,  and  the  final  three  levels  are  the  typical  three 
laser  levels  of  the  Cs  DPAL.  In  essence,  this  version  of  the  XPAL  simply  found  a  new 
way  to  populate  the  62P3/2  state  in  Cs.  A  temperature  of  435 °K  was  maintained  in 
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these  early  studies. 

In  2009,  Readle  et  al.  performed  a  similar  analysis  using  krypton  instead  of 
argon  [34],  This  time  the  Nd:YAG  pumped  dye  laser  was  tuned  to  841  nm  to  pump 
the  blue  satellite  peak  of  the  Cs-Kr  excimer.  The  Cs-Kr  was  contained  in  a  Pyrex 
cell  with  plane-parallel  windows,  having  a  length  of  10  cm  and  diameter  of  2.5  cm. 
The  cell  contained  500  Torr  of  Kr  and  100  Torr  of  C2H6.  The  emitted  laser  energy 
was  detected  by  calibrated  pyroelectric  detectors.  Three  temperatures  were  used  in 
this  experiment,  374  K ,  443  K  and  468  K.  The  spectral  breadth  of  the  excitation 
spectrum  blue  satellite  was  ^  2  nm.  This  is  smaller  than  the  Cs-Ar  blue  satellite,  but 
still  85  times  larger  than  the  typical  Cs  D2  line  in  the  DPAL  system.  It  was  found  that 
parasitic  lasing,  due  to  the  plane  parallel  windows,  contributed  to  a  decrease  in  laser 
performance  for  larger  alkali  densities.  Three  temperatures  and  their  corresponding 
thresholds  were  reported;  40  jiJ  at  444  K  (similar  to  the  previous  study),  85  ixJ  at 
458  K  and  140  /iJ  at  468  K.  It  was  expected  and  reported  that  both  the  threshold 
increased  and  the  slope  efficiency  decreased  with  increasing  alkali  number  densities. 

In  2010,  Readle  et  al.  changed  two  major  aspects  of  the  operation  of  the  XPAL 
system.  First  they  (which  caused  both  changes  to  occur  simultaneously)  removed 
the  hydrocarbon  from  the  equation.  This  was  done  because  of  the  degrading  effect 
hydrocarbons  has  on  lasing  mediums.  In  addition  to  a  very  narrow  spectral  pump 
bandwidth,  the  typical  DPAL  system  suffers  from  a  very  significant  degradation  in 
laser  performance  over  a  short  period  of  time  due  to  the  breaking  down  of  the  hy¬ 
drocarbon  in  the  alkali  cell.  The  breaking  down  of  the  hydrocarbon  creates  a  ”soot” 
in  the  lasing  medium  and  drastically  increases  the  loss.  By  removing  the  hydrocar¬ 
bon,  population  inversion  occurs  at  the  2Ps/2  state,  allowing  this  particular  XPAL  to 
operate  in  a  4-lcvel  system.  Test  cells  containing  Cs-Ar  and  Cs-Ar-Kr  mixtures  were 
explored.  Like  previous  studies,  500  Torr  of  alkali  were  used  in  both  test  cells.  In 
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the  case  of  the  Cs-Ar-Kr  cell,  equal  parts  of  each  rare  gas  were  used  until  500  Torr 
was  reached.  The  same  Nd:YAG  pumped  dye  laser  was  used  as  the  pump  source 
and  the  test  cell  was  10  cm  in  length.  The  breadth  of  the  excitation  spectrum  de¬ 
creased  slightly  from  the  Cs-Ar-C2H6  mixture  to  ^  3  nm.  It  was  also  noted  here  that 
the  spectral  width  narrows  slightly  at  higher  alkali  concentrations.  Surprisingly  the 
threshold  absorbed  pump  pulse  energy  at  464  K  was  ~  130  /i J,  which  was  very  close 
to  the  threshold  value  of  the  Cs-Kr-C2H6  laser  at  comparable  temperature.  As  in 
the  last  study,  three  temperatures  along  with  their  thresholds  were  reported  for  the 
Cs-Ar  4-level  system;  130  [iJ  at  464  K  (similar  to  the  468  K  case  in  previous  study), 
175  iiJ  at  475  K  and  250  / iJ  at  485  K.  It  was  noted  that  although  the  Cs-Ar-Kr 
mixture  had  the  same  spectral  blue  satellite  width  and  peak,  and  nearly  identical 
absorbed  pump  energy  threshold  conditions,  that  the  benefit  comes  by  reducing  the 
D2  laser  pump  energy  threshold  significantly  (as  much  as  28%). 

In  2011,  Palla  et  al.  presented  an  XPAL  model  and  theory  at  an  SPIE  confer¬ 
ence  [4],  A  BLAZE-V  model  was  modified  to  model  the  4-lcvcl  XPAL  system  for 
Cs-Ar  pairs.  The  model  was  validated  against  the  data  taken  in  the  previously  re¬ 
ported  studies  by  Reaclle  and  then  extended  to  predict  the  performance  for  a  Rb-Kr 
mixture.  The  model  suggested  that  the  Rb-Kr  variant  should  have  significantly  lower 
required  pump  intensities  than  the  Cs-Ar  cousin. 

In  2012,  Hewitt  et  al.  used  the  Cs-Ar  4-level  XPAL  to  explore  the  effect  photoion¬ 
ization  has  on  the  dynamics  of  the  laser  operations  when  using  pump  intensities  of 
up  to  3  [33].  Approximately  500  Torr  of  argon  was  used  in  a  10  cm  in  length 

borosilicate  cell.  Six  cell  temperatures  in  the  range  of  433  to  534  K  was  analyzed.  The 
following  non-linear  process  was  used  to  determine  how  significant  a  role  two-photon 
ionization  plays  in  XPAL  kinetics. 
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(4) 


(2) 

Cs(6p2P3/2 )  +  2hw  ►  Cs+(150)  +  e“  +  AA 

Where  AE  is  the  energy  beyond  that  required  for  two  836.7  nm  to  ionize  the  atom. 
For  a  474  K  case,  the  two-photon  ionization  cross-section  was  calculated  to  be  no 
larger  than  8  x  10”26  ^4  and  therefore  determined  that  photoionization  plays  only  a 
minor  role  in  the  laser  process  of  the  XPAL  even  at  these  high  pump  intensities. 

This  work  will  continue  the  investigation  of  the  XPAL  by  demonstrating  perfor¬ 
mance  of  the  Rb-Ar  and  Rb-Kr  4-level  system  in  a  heat  pipe.  A  heat  pipe  allows 
testing  the  performance  at  high  temperatures  and  high  pressures. 
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III.  The  Two-Photon  Absorption  (TPA)  Cross-Section  for 
the  5 2Si/2  — >  h2D5/2  Transition  in  Naturally  Occurring  Rb 


3.1  Introduction 

A  new  class  of  diode  pumped  lasers  based  on  excitation  of  the  D2  line  and  sub¬ 
sequent  lasing  on  the  Di  line  in  alkali  atoms  is  receiving  considerable  attention  [65]. 
The  Diode  Pumped  Alkali  Laser  (DPAL)  typically  requires:  (1)  high  buffer  gas  pres¬ 
sures,  >  1  atmosphere,  to  match  the  width  of  the  absorption  profile  to  the  spectral 
lineshape  of  the  diode  pump  source,  and  (2)  a  collision  partner  such  as  ethane  at 
pressures  of  several  hundred  Torr  to  relax  the  population  in  the  pumped  2P3/2  to  the 
upper  laser  2Pi/2.  The  alkali  atoms  provide  high  absorption  even  at  low  concentra¬ 
tion  and  a  single  atom  may  rapidly  cycle  under  intense  pump  conditions,  at  rates 
exceeding  1010  cyc^es ;  suggesting  high  power  operation  from  a  small  gain  volume. 

Several  additional  schemes  for  optically  pumping  alkali  metals  to  achieve  lasing 
in  the  infrared  and  blue  have  been  developed  over  the  past  40  years  [1,  5,  6,  7,  39, 
55,  59,  63].  More  recently,  stimulated  emission  on  the  blue  2P  — 2  S  transitions  in 
Rb  and  Cs  has  been  achieved  by  pumping  at  a  single  wavelength  in  the  red  via 
two-photon  absorption  [16].  The  52D3/2,  52D5/2,  and  72S'i/2  states  of  rubidium  and 
the  72D3/2,  72D5/2,  and  92S'i/2  states  of  cesium  were  populated  at  low  pressure  by 
two-photon  excitation  using  a  pulsed  dye  laser.  The  transition  probability  for  this 
two-photon  transition  from  the  initial  ground  state  | i),  or  52S'i/2,  to  the  final  state 
|/),  or  25 D5/2,  proceeds  through  a  virtual  intermediate  state  |u)  as  shown  in  Figure  7. 
The  pump  wavelength  is  12,851.749  cm-1,  which  is  35.204  cm-1  to  the  blue  of  the 
D2  line.  The  slope  efficiency  of  0.5%  in  Rb  is  considerably  higher  than  achieved 
in  the  prior  sequential  double  resonance  experiments  [63]  and  requires  only  a  single 
wavelength  pump  source.  Cascade  lasing  on  the  infrared  transition  followed  by  the 
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blue  transition  offers  no  quantum  defect  and  the  potential  for  low  heat  loads.  There 
would  also  be  no  need  for  a  spin-orbit  coupling  gas,  such  as  ethane,  which  has  caused 
soot  buildup  problems  in  traditional  DPAL  cells  at  high  temperature  [57].  We  refer 
to  this  blue  laser  as  the  Two-photon  Pumped  Alkali  Laser  (TPAL). 

The  blue  TPAL  is  also  an  example  of  a  rapid  non-linear  optical  process  that  can 
populate  higher  lying  states  in  the  alkali  atom.  This  is  of  interest  to  the  DPAL 
community  when  high  pump  intensities  are  used.  Emission  from  the  second  excited 
P  states,  and  from  more  than  30  lines  including  states  near  ionization  have  been 
observed  by  Jones  et  al.  [41].  If  the  ionization  rate  becomes  too  large,  the  situation 
could  impose  a  limit  to  intensity  scaling. 

crrr1 

2D5/2  - -  25,703.498 

2D3/2  25,700.536 


a- 778.1  nm 


2P3/2 

2Pl/2 


A  (detuning) 

12,851.749 

12,816.545 

12,578.95 


v=  778.1  nm 


2Sl/2  -  0 

Figure  7.  The  two-photon  absorption  cross-section  was  calculated  for  this  particular 
transition  in  Rb. 


Few  groups  have  modeled  or  measured  the  absorption  cross-section  for  the  52Si/2  — 
52D,5/2  transition  in  Rb.  Zapka  et  al.  used  Doppler-free,  frequency-modulation  (FM) 
spectroscopy  to  generate  the  absorption  spectra  in  Rb  vapor  in  1983  [67].  They 
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focused  a  60  mW  single-mode  cw  dye  laser  to  a  beam  waist  of  70  (irn  in  a  81  Rb  va¬ 
por  of  3.3  x  1014  cm-3  and  observed  an  absorption  of  2.5%.  A  TPA  cross-section 
(j 2(p  =  4  x  10~2°  ^  for  the  52Si/2(F  =  2)  — »  52D5/2(F  =  4)  hyperfine  line  in  8l Rb 
was  reported. 

Collins  used  Doppler  free,  differential  absorption  spectroscopy  (DAS)  [19]  with  an 
external-grating-cavity  cw  diode  laser  with  100  to  tune  through  the  52S'i/2(A  = 
2)  — »  52D5/2(F  =  0,1,  2,  3, 4)  transitions  of  85Rb  in  a  2.54  cm  long  Pyrex  cell 
containing  a  85 Rb  number  density  of  6.8  x  10”18  cm~3.  Five  Gaussians  were  fit  to 
the  weak  absorption  signals  (a  oc  10-6)  corresponding  to  the  five  hyperfine  lines  to 
determine  a  two-photon  absorption  cross-section  for  the  52>S'i/2  — »  52D5/2  transition 
of  cr 2<j>  =  1-2  x  10-18  was  two  orders  of  magnitude  larger  than  Zapka’s  report.  In 
1994,  a  theoretical  value  for  the  TPA  cross-section  was  calculated  (0.57  x  10-18  A%) 
using  reduced  absorption  rates  [49].  Most  recently,  TPA  cross-sections  were  measured 
using  low  powers  (1  mW )  and  low  densities  in  a  hollow-core  photonic-band-gap  fiber 
(PBGFs)  filled  with  Rb  vapor  [43].  In  this  research  PBGFs  were  demonstrated  to 
enhance  the  light-matter  interaction  for  nltralow  powers,  mitigating  the  need  for 
focused  beams  and  high  powers.  A  1%  absorption  was  determined  by  detecting 
30  nW  of  output  blue  flourescence  at  420.3  nm  with  a  1  mW  diode  pump  laser.  With 
this  method  a  direct  absorption  value  was  calculated  to  be  5  x  10~20  for  the 
52S1/2  — >•  52D5/2  in  S5Rb. 

There  has  not  been  any  consistency  in  the  values  across  the  four  previous  groups. 
Therefore,  a  more  comprehensive  study  is  required  to  better  understand  the  true 
value.  However,  generating  the  needed  high  intensity  is  problematic.  Bamford  et  al. 
reported  a  technique  using  pulsed  lasers  for  determining  the  two-photon  absorption 
cross-section  in  atomic  oxygen  [20].  Peak  intensities  of  1.8  were  achieved,  but 
a  careful  determination  of  both  the  spatial  and  temporal  profiles  was  required  to 
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accurately  extract  a  cross-section.  Modeling  the  temporal  profile  is  also  cumbersome 
since  the  run-to-run  variation  of  the  intensity  profile  is  high.  Therefore  our  research 
chose  to  use  a  focused  cw  approach  that  varied  the  Rb  number  density  from  1.61  — 
8.52  x  1015cm-3,  ranged  pump  powers  from  0.32  —  0.99  W  (actual  power  measured 
at  test  cell)  and  used  two  focal  length  lenses,  15.24  cm  and  20.32  cm,  to  generate  87 
test  configurations.  It  was  necessary  to  develop  a  model  of  the  spatial  profile  for  the 
intensity  to  properly  interpret  the  data.  This  model  consists  of  a  linear  combination 
of  the  two-photon  Beer’s  Law  variant  and  a  Gaussian  propagation  term  to  account 
for  the  focused  beam.  The  statistical  average  of  the  87  TPA  cross-sections  for  the 
52S'i/2  — >  52D5/2  transition  in  naturally  occurring  Rb  will  provide  a  confident  value 
to  be  used  to  support  work  in  examining  the  pump  rate  for  a  TPAL  system. 

3.2  Experimental 

Figure  8  illustrates  the  major  elements  of  the  experimental  set-up.  Although  a 
2.5  W  Sirah  Matisse  cw  Tksapphire  ring  laser  was  used  as  the  pump  source  the  max 
output  power  did  not  exceed  1.2  W.  The  Matisse  stabilized  the  output  power  as  it 
scanned  in  frequency  by  10  —  15  GHz  using  a  Pound-Drever-Hall  waveform  technique. 
Typical  scan  rates  were  set  at  50  About  1%  of  the  power  was  sent  to  a  HighFi- 

nesse  WS/7  High  Precision  wavemeter  to  determine  the  frequency  to  within  60  MHz 
accuracy  on  a  300  THz  scale.  The  beam  then  passed  through  a  mechanical  chopper 
at  2  kHz  frequency  and  then  to  a  high  reflector  (Ml).  Another  1%  of  the  beam 
power  is  used  to  monitor  the  incident  power  and  transmission  through  a  300  MHz 
etalon,  while  the  remaining  energy  is  sent  through  a  lens  to  focus  the  beam  into  the 
alkali  test  cell.  As  depicted  in  Figure  8,  two  Thor  Lab  PDA55  photodiodes  were  used 
to  monitor  the  incident  and  transmitted  power.  After  exiting  the  test  cell  the  beam 
passed  through  a  nuetral  density  filter  wheel,  OD  =  0. 1-5.0,  to  prevent  detector  satu- 
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ration.  The  signals  from  the  two  detectors,  D]  and  D2,  and  the  etalon  each  were  sent 
through  their  own  Stanford  Research  System  SR850  lock-in  preamplifier  triggered  by 
the  chopper  wheel.  Integration  times  of  up  to  1  s  were  used.  The  raw  modulated 
signals  were  maintained  and  packaged  into  LabView  files  via  a  National  Instruments 
USB-6212  BNC  multifunction  data  acquisition  (DAQ)  module  capable  of  400 

The  alkali  pyrex  cell  was  15  cm  in  length  and  was  fixed  to  a  gas  handling  and 
vacuum  system  with  an  ultimate  vacuum  of  29  mTorr.  The  alkali  cell  was  enclosed 
inside  a  metal  oven.  The  temperature  of  the  oven  was  kept  at  about  10  degrees  hotter 
than  the  finger  containing  the  solid  alkali  to  prevent  condensing  on  the  windows. 
Three  temperature  zones  were  controlled  and  monitored  by  an  Omega  Engineering 
temperature  controller.  The  first  two  zones  are  along  the  length  of  the  test  cell 
(opposite  ends)  and  the  third  on  the  finger  protruding  under  the  center  of  the  cell. 


Figure  8.  Two-photon  cw  laser  absorption  apparatus. 


A  300  M Hz  Coherent  Laser  Spectrum  Analyzer  was  used  to  assist  in  determining 
the  laser  frequency.  We  self  calibrated  the  free  spectral  range  (FSR)  at  302.5  MHz 
using  the  hyperfine  structure  of  the  ground  state  in  Rb. 
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The  power  just  inside  the  test  cell  was  measured  by  placing  a  Coherent  LabMax- 
TO  power  meter  via  the  LM-3  HTD  power  sensor  in  the  optical  beam  train  where 
the  alkali  test  cell  would  normally  be.  An  exact  replica  of  the  1”  Pyrex  glass  window 
was  added  into  this  beam  train  to  simulate  the  the  power  entering  the  absorption 
medium.  Transmission  issues  were  encountered  in  the  test  cell  when  operating  at 
temperatures  from  280  °C  to  320  °C.  Liquid  droplets  of  alkali  were  observed  raining 
down  from  the  neck  to  the  test  cell.  This  disturbance  started  to  become  noticeable 
at  280  °C  and  quite  dominant  at  temperatures  above  300  °C.  Therefore,  all  of  the 
TPA  scans  used  in  this  research  were  between  230  and  280  °C.  These  temperatures 
were  measured  by  monitoring  the  finger  temperature  of  the  alkali  with  an  RTD  and 
the  Omega  Temperature  Controller.  Since  the  vapor  pressure  is  very  sensitive  to  the 
temperature  of  the  alkali  finger,  care  was  given  to  understanding  the  concentration 
of  the  alkali  vapor  inside  the  test  cell. 

Absorption  in  the  wings  of  the  optically  thick  D2  transition  was  used  to  monitor 
the  Rb  concentration.  Figure  9  illustrates  an  example  spectrum.  The  data  is  rep¬ 
resented  in  blue,  while  the  red  line  is  the  model.  The  first  peak  contains  the  three 
transitions  of  5 2S1/2(F"  =  1)  — ►  5 2P3/2(F'  =  0),  5 2S1/2{F"  =  1)  — ►  5 2P3/2(F'  =  1) 
and  52Si/2(F"  =  1)  — »  52P3/2(F'  =  2)  of  87Rb.  At  85  °C,  the  peak  absorbance  is 
very  large  ( A  >  22).  A  fit  of  the  full  Voigt  broadened  hyperfine  spectra  is  used  to 
determine  the  Rb  concentration. 

A  =  ~ln  ({)  =  G(v'  Vr"^F,)  (fy)  (sf)  AfiNL  (5) 

where  G(u,  vF"^f')  —  l'F"^F,)SF"F,fBfiso  is  the  hyperfine  lineshape 

F'  F"  iso 

function  that  accounts  for  the  hyperfine  transitions  F”  — >  F',  are  the  line 

strengths,  /#  is  the  Boltzmann  factor,  fjso  is  the  isotopic  factor,  gj  is  the  degener¬ 
acy,  and  L  is  the  length  of  the  test  cell.  Between  the  two  isotopes,  there  are  12  allowed 
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hyperfine  transitions.  The  model  constrains  the  known  spacing  and  the  relative  am¬ 
plitudes  of  each  hyperfine  component.  These  amplitudes  were  calculated  from  the 
line  strengths,  Sf"f >,  the  Boltzmann  factor,  and  the  isotopic  fraction,  fiso. 

Seven  parameters  were  allowed  to  float;  the  total  Rb  concentration,  the  frequency  of 
the  first  line,  the  Doppler  width,  and  4  parameters  for  the  cubic  baseline.  In  Figures  9 
the  temperature  was  measured  by  the  RTD  as  85  °C,  with  a  pressure  that  measured 
at  34  uTorr. 


Figure  9.  The  entire  D2  transition  is  centered  at  approximately  384.230  THz.  Although 
the  intensity  was  kept  below  Isati  the  high  temperature  caused  absorbances  to  exceed 
the  dynamic  range  of  the  SR850  lock-in.  The  pressure  and  temperature  for  this  run 
were  measured  at  39.1  mTorr  and  85  °C  respectively.  The  blue  line  refers  to  the  actual 
data  collected  and  the  red  is  the  Voigt  model  fitted  to  the  transition. 


The  spectroscopically  determined  Rb  vapor  pressures  are  compared  with  a  known 
vapor  pressure  curve  (dotted  line)  [62]  in  Figure  10.  A  curve,  using  the  same  pa¬ 
rameters  from  the  established  vapor  pressure  curve,  was  shifted  in  temperature  to 
overlay  our  data.  This  shift  in  temperature  (shown  in  Figure  10)  amounted  to  about 
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14  ±  0.5  °C  cooler  than  measured  at  the  walls  of  our  alkali  test  cell.  The  0.5  °C 
uncertainty  in  temperature  corresponds  to  an  11%  change  in  alkali  concentration. 


60  70  80  90  100  110  120  130  140  150  160 

Temperature  (degrees  Celcius) 


Figure  10.  The  vapor  pressure  calibration  curve.  The  temperature  inside  the  alkali 
test  cell  was  approximately  14  °C  cooler  than  measured  on  the  wall  of  the  cell. 

The  beam  quality,  or  M2,  of  the  laser  beam  was  measured  using  a  ModeMaster 
by  Coherent.  A  range  of  powers  were  used  from  1090  mW  to  1320  rnW  (measured  at 
the  pump  laser  aperture),  to  understand  the  power  dependence  on  beam  quality.  A 
filter  wheel  would  have  affected  the  quality  of  the  beam  and  was  not  used.  Instead, 
the  pump  beam  to  the  Tksapphire  was  adjusted.  Referring  to  Figure  8,  M2  was 
measured  at  the  location  of  the  alkali  test  cell.  It  was  observed  that  over  this  power 
range  the  beam  quality  stayed  fairly  consistent,  from  1.05  —  1.09  ±  0.03  .  We  chose 
M2  =  1.09  to  be  conservative  while  calculating  the  final  number  for  a^- 

A  ModeMaster  was  again  placed  in  the  beam  train  where  the  test  cell  would 
normally  be,  minus  the  lens.  At  1200  mW  the  ModeMaster  measured  both  the 
external  beam  waist,  a jext  =  0.547  mm,  and  the  distance  between  the  external  waist 
to  the  ModeMaster  position,  z\  =  60  cm.  These  quantities  allowed  us  to  calculate 
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the  beam  waist  at  the  focus; 
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i--Y  + 
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V/AA/V 


(6) 


Two  lenses  were  used  (15.24  cm  and  20.32  cm)  to  investigate  the  veracity  of 
this  work’s  ability  to  extract  a  consistent  TPA  cross-section  number,  regardless  of 
a  change  in  intensity  due  to  the  beam  waist.  These  lenses  were  inserted  into  the 
beam  train  such  that  the  focus  landed  at  the  center  of  the  alkali  test  cell.  This  was 
done  for  symmetry  purposes  that  will  become  clear  during  the  development  of  the 
“Governing  Equation”  which  modeled  a  focused  Gaussian  beam  propagating  through 
a  two-photon  absorption  medium. 


3.3  Results 

The  cw  laser  absorption  spectrum  for  the  52S'i/2  — >  52D5/2  two-photon  transition 
for  an  intensity  at  the  waist  of  /  =  3.1  YiT  ,  Rb  concentration  of  N  =  4.07  x  1015  ator"ls 
and  the  test  cell  path  length  of  L  =  15  cm  is  shown  in  Figure  11. 

For  a  collimated  beam,  the  transmitted  intensity  decays  with  propagation  dis¬ 
tance,  z,  as: 


l  -  ~^N  <7> 

jm~l  =  a*tNI° L  =  A  (8> 

where  a2(j,  is  the  two-photon  absorption  cross-section  (q^)  and  L  is  the  propagation 
distance  in  the  context  of  the  length  of  the  test  cell.  For  a  focused  beam,  the  Jq 


factor  in  the  absorbance  might  be  best  interpreted  as  a  weighted  average  of  the 
intensity.  A  more  complete  development  of  the  spatial  dependence  of  the  intensity 


Figure  11.  A  typical  spectra  collected  in  this  experiment.  Peak  (1)  includes  hyperline 
transitions  of  two  isotopes.  Peak  (2)  was  used  throughout  this  study  to  measure  the 
TPA  cross-section  due  to  it’s  independence  from  other  isotopes  and  relative  strength 
compared  to  Peak  (3).  The  beam  waist  was  calculated  to  be  69.7  /i m .  The  vertical  lines 
refer  to  the  theoretical  positions  of  the  hyper-fine  splittings  of  the  ground  state  of  each 
Rb  isotope.  The  outer  two  lines  are  Sl  Rb  and  the  inner  two  are  85 RB. 

will  be  developed  below.  Integrating  the  intensity  over  the  detector  area  yields  a 
signal  proportional  to  the  power,  P ,  as  presented  in  Figure  11.  A  peak  absorbance  of 
A  =  0.19  is  observed.  The  baseline  noise  is  ~  A  =  0.002,  with  a  longer  term  baseline 
variance  of  up  to  0.01. 

The  hyperfine  splitting  of  the  ground  52S'i/2  state  is  ~  3.036  and  6.834  GHz  for 
the  85 Rb  and  8l  Rb  isotopes  respectively  [62],  The  hyperfine  splittings  for  the  upper 
52_D,5/2  are  smaller  than  the  Doppler  width,  2.9  —  9.4  and  15.9  —  28.8  MHz,  and 
unresolved  (refer  to  Figure  12)  [3].  The  less  abundant  (27.83%)  8l  Rb  isotope  with 
the  larger  ground  state  hyperfine  splitting  has  7,  F  —  2, 1,  0  allowed  transitions  in 
two  groups.  The  more  abundant  (72.17%)  85 Rb  isotope  with  the  smaller  hyperfine 
splitting  has  10  allowed  transitions,  with  one  of  the  two  groups  convolved  with  the 
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low  frequency  component  for  8' Rb.  The  isotopic  shift  is  small,  «  78  MHz  [62],  The 
frequency  differences  observed  in  Figure  11  are  reduced  by  a  factor  of  two,  due  to  the 
two-photon  transition. 


F=1 

F=2 

F=3 


F=4 


85Rb 


87Rb 


Figure  12.  The  hyperflne  splittings  for  both  isotopes  [3]. 


We  chose  to  report  values  from  the  52Si/2(F”  =  2)  — *  52D5/2(F'  =  0,1,  2,  3, 4) 
transition  in  85 Rb  due  to  its  signal  strength  and  independence  from  the  other  isotope. 
This  transition  is  labeled  as  Peak  #2  in  Figure  11.  With  the  absorbance  of  this 
peak  equal  to  be  0.115,  we  calculate  the  effective  “one-photon”  cross-section  to  be 
(<720  Iq)  ~  1.88  x  10-18  cm2  from  Equation  8.  The  difficulty,  and  the  focus  of  this 
research,  is  determining  with  accuracy  the  intensity  through  the  alkali  cell  such  that 
the  TPA  cross-section  can  be  extracted.  To  do  this,  both  a  spatial  model  of  the 
intensity  propagating  through  the  alkali  cell  and  a  lineshape  model  are  required. 

It  is  important  to  understand  the  TPA  cross-section  number  in  the  context  of 
the  lineshape,  because  the  absorption  cross-section  is  directly  proportional  to  the 
lineshape.  The  lineshape  model  in  this  experiment  consists  of  a  determination  of 
the  Boltzmann  and  isotropic  factors,  the  two-photon  signal  strengths,  the  frequency 
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Table  2.  The  calculated  signal  strengths  for  85  rubidium 
where  I  =  |. 


F’ 

0 

1 

2 

3 

4 

5 

2 

1/36 

7/90 

1/9 

7/60 

1/12 

0 

3 

0 

1/180 

1/36 

7/90 

1/6 

11/36 

spacing  and  a  lineshape  function.  For  this  transition,  two  Boltzmann  factors  were 
calculated,  one  for  each  of  the  two  hyperhne  lines  of  the  ground  state  ( F ").  These 
two  values  were  calculated  to  be  /b(f"= 2)  =  0.417  and  fB(F"= 3)  =  0.583.  The  isotopic 
ratios  are,  flso  ssRb  =  0.7217  and  flS0  s7Rb  =  0.2784.  The  17  relative  signal  strengths, 
Sf"f',  for  a  two-photon  transition  were  calculated  using  the  following  equation  [24], 


(2  F”  +  1)(2F'  +  1) 


Sf'\f>  — 


(2/  +  1) 


(9) 


Tables  2  and  3  list  the  17  signal  strengths  calculated  from  Equation  9.  Finally,  the 
frequency  separation  was  determined  from  data  provided  by  both  Steck  and  Olsen 
[3,  62],  Figure  11  shows  the  theoretical  line  spacings,  however  the  F  spacings  are 
unresolvable  on  that  scale.  We  chose  an  area  normalized  Voigt  function  to  model 
the  absorption  lineshapes  due  to  the  need  to  model  the  Doppler  and  homogeneous 
linewidth.  We  discovered  that  there  was  no  preferred  function  used  to  model  a  two- 
photon  transition  from  previous  studies  [2,  10,  19,  21,  43,  45].  Although  we  used  a 
Voigt  function,  the  Doppler  width  is  approximately  100  times  larger  than  the  homo¬ 
geneous  width,  so  the  Gaussian  lineshape  is  dominant. 

Before  applying  a  lineshape  model  to  the  data  in  Figure  11,  a  non-linear  transfor¬ 
mation  of  the  y-axis  from  absorbance  units  to  cross-section  units  is  performed.  A  full 
interpretation  of  the  spectra  to  determine  a  best  estimate  of  the  two-photon  cross- 
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Table  3.  The  calculated  signal  strengths  for  87  rubidium 
where  I  =  |. 


F’ 

1 

2 

3 

4 

1 

9/80 

7/48 

7/60 

0 

2 

1/80 

1/16 

7/40 

3/8 

section  requires  an  analysis  of  the  focused  beam  propagation.  Figure  13  illustrates 
the  intensity  propagation  for  a  focused  beam. 

A 


l 


Figure  13.  The  laser  beam  passes  through  a  focus  as  it  travels  through  the  Rb  Vapor. 
An  accurate  calculation  of  the  TPA  cross-section  requires  a  model  of  I(z). 


In  differential  form,  this  propagation  is  modeled  by  the  following  equation, 


dl 

dz 


-aNI2  +  F(z)I 


where  the  Gaussian  beam  propagation  term  is: 


(10) 


F(z) 


2  (  r2  \duj 

l o(z)  \o;2(z)  /  dz 


(a) 
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where, 


uj(z)  =  ul 


1  +  a  (z  -  — ) 


and  the  parameters  are  define  as: 


(12) 


M2  A 
a  = - = 

™ o 

tao  =  beam  waist  =  69.7  —  94.4  ym 
M 2  =  beam  quality  factor  =  1.09 


(13) 

(14) 

(15) 


The  solution  to  Equation  10  follows  from  Bernoulli’s  approach,  such  that: 


e[/i(L)-/i(0)] 

UL)  =  S  +  ff2#lv/„[/2(i)-/2(o)]  <?-'.«> 


(16) 


where 


and, 


h(L) 


2o?(L~H)  (  r2  \  dL 

1  +  a\L  -  Hf2  [1  +  a2(L  -  H)2}  ) 


(17) 


I2{L)  =  /  eh{L)dL 


(18) 


Unfortunately,  the  integral  1-2  is  not  analytic  and  a  numerical  approach  is  required. 
We  choose  to  cast  the  solution  in  the  form  of  Equation  8: 


m  = 


1  +  <T2(j)Iov(aJ  r)NL 
where  the  longitudinally  averaged  intensity  is  developed  from  the  factor: 


(19) 
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V(<*,r)  =  ~j~ 
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r2 

~Z72 


'6\i+^(*-§)2Jdz 
1  +  a2(z  —  j)2 
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a2L2  z? 

(!+— )e 


6  V  1+^^ 


(20) 


A  more  complete  development  of  Equations  16—20  is  found  in  Appendix  A.  The 
average  intensity  for  two  Gaussian  beams  is  shown  as  a  function  of  radial  position  in 
Figure  14.  As  the  beam  becomes  more  collimated,  a  — >  0,  the  solution  becomes  less 
sensitive  to  radial  position  and  the  solution  of  Equation  8  is  recovered. 


Figure  14.  The  intensity  varies  as  a  function  of  rj(r)  in  the  radial  direction.  As  the 
beam  waists  decreases,  or  a  — >  oo,  the  rj  function  plays  a  more  significant  role  in  the 
intensity  equation. 


Figure  15  illustrates  the  longitudinal  intensity  distribution  on  axis  for  several  TPA 
cross-sections.  The  first  is  an  unattenuated  case,  or  a  =  0,  to  set  the  baseline  and 
contrast  the  effects  of  TPA.  The  next  three  are  in  powers  of  ten,  cr  =  10-21, 10~20  and 
10~19.  The  location  of  the  peak  intensity  shifts  to  the  left  with  increasing  TPA  cross- 
section.  This  suggests  that  the  TPA  term  in  Equation  10  dominates  the  focusing 
term  for  certain  values  of 
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cr  =  0 


Figure  15.  Using  Equation  19,  this  plot  depicts  four  intensity  profiles  corresponding 
to  four  different  TPA  cross-section  values.  Since  the  first  value  is  zero,  we  find  that 
the  equation  returns  a  symmetric  Guassian  as  expected.  Please  note  that  as  the  cross- 
section  increases,  that  the  intensity  at  the  end  of  the  test  cell  (L  =  15  cm)  decreases 
and  the  peak  of  the  intensity  profile  increasingly  shifts  to  the  left  of  the  focus. 


Figure  16  illustrates  three  methods  for  modeling  the  absorbance  value.  The  y-axis 
is  the  typical  Beers  Law  relation  for  the  absorbance,  where  three  methods  have 

been  fit  to  this  form.  The  three  methods  represented  on  this  plot  are  the  One-Photon 
Beers  Law  (Equation  21),  the  two-photon  Beers  Law  (Equation  22)  and  the  two- 
photon  Beers  law  with  a  focused  Gaussian  beam  (Equation  23). 


—In^r  =  a2(j>IQNL 
h 

(21) 

In^r  =  ln(  1  +  a2(t>I0NL) 

(22) 

=  ln(  1  +  a2rj>r](a,r)I0NL) 

(23) 

Each  method  is  plotted  against  the  two-photon  absorbance,  A2<j,  =  cj^IqNL ,  along 
the  x-axis.  Comparing  the  One-Photon  Beers  Law  (1)  to  the  two-photon  Beers  Law 
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(2),  we  discover  that  as  the  absorbance  increases,  the  One-Photon  Beers  Law  overesti¬ 
mates  the  value  for  the  TPA  cross-section.  Finally,  we  see  in  Figure  16  how  the  model 
for  the  two-photon  Beers  Law  with  a  focused  Gaussian  beam  (3)  compares  to  the  two- 
photon  Beers  Law  (2)  for  two  beam  waists  (a  =  0.303  and  0.555).  Interestingly,  for 
the  smaller  beam  waist,  both  the  One-  and  two-photon  Beers  Laws  underestimate  the 
value  of  the  TPA  cross-section  for  low  values  of  absorbance  and  then  overestimates 
by  an  increasing  amount  for  larger  values  of  absorbance.  We  see  that  for  a  =  0.555, 
the  two-photon  Beers  Law  is  a  good  approximation  for  G-2<t>I0NL  <  0.23.  Figure  16 
also  shows  that  for  larger  beam  waists  {i.e.  a  =  0.303)  and  larger  absorbance  values, 
there  is  greater  divergence  between  the  two-photon  Beers  law  model  and  the  numerical 
model  developed  in  this  work.  This  divergence  points  to  the  need  to  properly  model 
the  pump  beam  going  through  a  focus  to  more  accurately  predict  the  absorbance 
value  for  a  given  TPA  cross-section. 

A  key  step  to  developing  a  function  that  transforms  absorbance  data  from  Fig¬ 
ure  11  to  TPA  cross-sectional  units,  is  to  first  perform  a  radial  integral  to  the  intensity 
model  shown  in  Figure  15  and  Equation  19.  This  radial  integral,  solved  numerically, 
illustrates  the  power  transmitted  through  the  alkali  cell.  This  integral  was  solved 
multiple  times  for  a  given  range  of  TPA  cross-section  values.  Next,  we  calculate  a 
two-photon  absorbance  value  from  the  computed  range  of  powers.  Plotting  these  cal¬ 
culated  absorbance  values  (i.e.  the  predicted  absorbance)  against  the  range  of  TPA 
cross-sections  yields  the  curve  in  Figure  ??.  A  4th-order  polynomial  is  then  fitted 
over  this  curve  to  extract  the  function  that  will  convert  (or  transform)  the  absorbance 
values  along  the  y-axis  in  Figure  11  to  TPA  cross-section  values.  This  transformation 
was  performed  on  all  87  test  configurations. 

After  transforming  the  data  in  Figure  11,  the  y-axis  now  represents  the  absolute 
two-photon  absorption  cross-section  which  takes  into  account  the  hyperfme  structure 
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Figure  16.  Three  models  for  the  intensity  plotted  against  the  two-photon  absorbance. 
(1)  Models  the  intensity  using  the  one-photon  Beers  Law  standard  model.  This  model 
was  baselined  by  forcing  the  other  two  models  to  equal  —In  j  This  is  why  this  model 

appears  as  a  perfect  45°  line.  (2)  This  is  the  two-photon  Beers  law  model.  (3)  and 
(4)  are  both  the  complete  numerical  method,  taking  into  account  the  focused  Gaussian 
beam  and  integrating  over  all  of  the  radial  length.  Two  curves  appear  because  two 
different  beam  waists  were  used  in  this  experiment.  All  4  curves  were  generated  from 
Equations  21  to  23. 


of  the  transitions,  cr^F .  Using  the  following  definition  for  the  one-photon  absolute 
absorption  cross-section, 


&i(t>(y)  —  ^Afig(u1vF"^F')  (24) 

the  one-photon  absolute  cross-section  that  accounts  for  the  hyperhne  structure,  cr^F, 
can  be  written  using  the  hyperhne  lineshape  function  defined  in  Equation  5, 

aif(u)  =  ^rAfiG(u’uF"^F')  (25) 
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Therefore,  applying  Equation  73  to  Equation  25,  the  hyperfine  cross-section  in  terms 
of  the  cross-section  is  now, 


G(lS,  l/F"_yF>) 

g(v, 


At  line  center,  u0l  Equation  26  can  be  further  simplified  to 


(26) 


yy  (md)  =  ^(t'o)  VVV  Sf"  f1  f b(.F")  fiSO  (27) 

F'  F"  iso 

Equation  74  is  directly  applied  to  the  two-photon  case  by  making  the  follow¬ 
ing  substitutions,  cr^F(i/0)  — *  (vo)  and  cti^uq)  — >  020(^0) •  A  Voigt  model 

(based  on  Equation  74)  was  developed  to  fit  over  the  data  collected  in  this  experi¬ 
ment  as  shown  in  Figure  17.  After  the  y-axis  has  been  transformed  to  TPA  units, 
the  amplitude  of  the  Voigt  function  in  the  model  (at  line  center)  yields  the  hyperfine 
cross-section,  (z'o)-  The  hyperfine  independent  TPA  cross-section  at  line  center, 
020 (z'o),  is  being  defined  as 

^  =  (|)  64 <28) 

Where  Xp  is  the  pump  wavelength  (\p  =  778.104  rim) ,  h  is  Planck’s  constant  in 
[J  —  s],  A  is  the  detuning  of  the  virtual  state  from  the  52Ps/2  state  in  [1/s]  (as 
depicted  in  Figure  7)  and  is  an  effective  decay  rate.  Selection  rules  dictate  this 
transition  (D-S)  is  forbidden  for  a  single  photon,  so  this  decay  rate  term  does  not 
describe  a  photon  from  the  D-level  decaying  down  to  the  ground  state.  It  is  rather 
a  convenient  collection  of  terms  that  include  the  two  dipole  moments  (/i 32  and  /i2i) 
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that  significantly  contribute  to  the  transition, 


20  64  7 tavI  2  2 

7/!  = 


(29) 


With  these  terms  defined,  the  y-axis  transformed  and  the  Voigt  model  developed, 
a  final  analysis  is  now  done  of  the  data.  An  example  of  this  final  fit  is  shown  in 
Figure  17. 
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Figure  17.  An  example  of  the  Voigt  model  fit  over  the  TPA  spectra.  The  amplitude  of 
the  Voight  was  extracted  to  determine  the  TPA  cross-section. 


This  final  fit  was  done  on  all  87  test  configurations,  allowing  87  TPA  cross-sections 
to  be  extracted.  These  TPA  cross-section  values  are  displayed  in  Figure  18  and 
compared  against  two  historical  values,  as  well  as  the  value  calculated  from  Equa¬ 
tion  74.  This  work  reports  an  average  experimental  value  for  the  hyperfine  broadened 
TPA  cross-section  for  the  52Si/2  — >  52D5/2  transition  in  natural  rubidium,  weighted 
against  the  error  of  each  test  case,  to  be  6.75  x  10-21  We  also  report  a  standard 
deviation  of  3.61  x  10~21  The  theoretical  value,  which  accounts  for  the  dipole 
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transition  moments  (Equation  74)  is  2.50  x  10  21 
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Figure  18.  All  87  TPA  cross-section  values  compared  against  historical  values.  Zapka 
et  al.  [67]  and  Saha  et  al.  [43]  are  within  a  factor  of  10  higher  in  cross-section,  while 
Marinescu  et  al.  (ct20  =  5.7  x  10-17e(^)  [49]  and  Collins  et  al  (020  =  1.2  x  I0~ls^r) 

[19]  are  off  the  chart.  Also  displayed  is  the  weighted  average  computed  in  this  work, 
along  with  the  value  computed  by  the  dipole  moment  dependent  equation  derived  from 
theory  (Equation  28). 


3.4  Discussion 

There  are  two  main  areas  that  contribute  to  significant  error  in  our  number.  The 
first  is  the  measurement  of  the  alkali  vapor  pressure  curve.  It  was  determined  that 
the  true  alkali  concentration  was  14  ±  0.5  °C  cooler  than  what  was  measured  by  the 
RTD  on  the  finger.  We  also  chose  to  simplify  the  result  by  assuming  no  significant 
temperature  gradients  existed  inside  the  alkali  test  cell.  We  felt  that  this  was  a 
reasonable  approximation  due  to  the  small  volume  of  the  test  cell.  A  second  significant 
factor  in  the  error  of  the  final  TPA  cross-section  value  was  the  determination  of  the 
spot  size.  The  primary  method  of  measuring  the  beam  size  was  to  use  a  ModeMaster 
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to  measure  the  beam  size  at  the  lens  position,  u>ext,  and  the  beam  quality,  M2.  The 
beam  size  was  measured  to  an  accuracy  of  ±  0.15  mm  and  the  beam  quality  to  within 
±  0.4. 

The  numbers  calculated  in  this  experiment  and  the  historical  values  must  be 
compared  only  in  the  context  of  their  respective  lineshape  functions.  All  previous 
experiments  utilized  sub-Doppler  techniques  to  extract  a  TPA  cross-section  value. 
First  of  all,  the  values  determined  by  Collins  and  Marinescu  are  greater  than  8000  x 
and  175  x  our  reported  value,  respectively.  Although  it  is  difficult  to  know  why 
these  two  groups  differ  so  greatly,  it  is  probable  that  the  source  of  the  error  lies 
in  the  assumptions  made  for  the  intensity  profile  and  the  ability  to  know  the  alkali 
concentration.  Collins  et  al.  also  did  not  do  a  direct  absorption  measurement  of 
the  transition.  Curves  generated  from  his  work  were  collected  from  a  side-viewing 
PMT,  so  an  incomplete  picture  may  have  resulted  from  not  being  able  to  compare  the 
initial  to  transmitted  intensity  properly.  Marinescu’s  work  centers  the  theory  around 
reduced  absorption  rates.  A  completely  divergent  method  from  the  rest  of  the  studies. 
The  values  reported  by  Zapka  and  Saha  are  only  6x  and  7x  our  reported  value, 
respectively.  We  believe  that  the  difference  between  the  line  widths,  and  ultimately 
the  lincshapes,  are  what  accounts  for  the  dissimilarity  in  TPA  cross-sections.  To  prove 
this  theory,  we  first  computed  the  natural  lineshape  based  on  the  widths  provided  in 
their  studies.  Dividing  their  reported  value  for  the  cross-section  by  these  lineshapes 
and  then  multiplying  by  our  Doppler  lineshape  of  similar  conditions,  we  calculate 
numbers  that  nearly  match  our  reported  weighted  average  value.  This  is  displayed  in 
both  Figure  18  and  in  Table  4.  In  an  effort  to  further  compare  our  value  with  those 
reported  in  the  literature,  we  then  applied  our  intensity  model  to  the  conditions 
under  which  Zapka  et  al.  observed  their  spectra.  We  chose  to  use  Zapka’s  reported 
conditions  primarily  because  they  were  the  most  similar  to  our  testing  conditions. 
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Table  4.  Comparing  the  previously  reported  values  to  our 
weighted  average  value  for  the  Two-Photon  Absorption 

(TPA)  cross-sections  measured  for  the  52<S'i/2  - >  52D.5/2 

Transition  in  Rb. 


Doppler-Free 

(  cm 4  \ 

aW  {~\ V) 

Reference 

Year 

Type 

Isotope 

Doppler-Broadened 

(  cm4  \ 

\lv) 

4xlO-20 

Zapka  et  al  [67] 

1983 

E 

87Rb 

1.30X10"21 

1.2xl0~18 

Collins  et  al  [19] 

1993 

E 

85Rb 

5.58  xlO"20 

0.57xl0~18 

Marinescu  et  al  [49] 

1993 

T 

85Rb 

N/A 

6.5xlO“20 

Saha  et  al  [43] 

2011 

T 

85Rb 

3.51xl0-21 

5  xlO~20 

Saha  et  al  [43] 

2011 

E 

85Rb 

2.70xl0-21 

N/A 

This  Work 

2013 

E 

85Rb 

6.78xl0-21 

Both  experiments  employed  a  cw  laser  and  a  lens  that  resulted  in  a  beam  waist  of 
~  70  fim.  Due  to  the  similarity  of  the  set-up  presented  in  his  work,  we  decided  to 
apply  our  numerical  model  to  the  conditions  Zapka  reported.  Our  model  calculated 
a  TPA  cross-section  of  3.5  x  10~20  g^,  which  was  within  12%  of  his  reported  number 
of  4.0  x  10^20  .  This  result  not  only  proved  greater  confidence  in  our  model, 

but  also  in  the  accuracy  of  Zapka’s  reported  result.  Next,  we  applied  the  theoretical 
equation  presented  in  Equation  74  to  calculate  a  TPA  cross-section  using  the  50  M Hz 
Lorentzian  line  width  reported  in  the  Saha  et  al.  work.  That  number  came  to  within  a 
factor  of  2  to  Saha’s  reported  value.  Although  Saha  et  al.  was  the  first  to  use  a  PBGF 
to  detect  this  weak  transition  using  ultra-low  powers,  many  similarities  existed  such 
as  focusing  the  beam  down  to  a  small  beam  waist  (~  6 urn,  10  x  smaller  than  Zapka 
and  ours)  and  collecting  an  input  and  output  intensity  signal  with  a  photodiode.  For 
a  complete  list  of  the  historical  values,  as  reported  in  the  literature,  please  refer  to 
Table  4. 

Having  the  TPA  cross-section  known  for  this  particular  transition  allows  for  an 
approximation  to  be  made  about  the  pump  intensity  threshold  condition  for  a  poten¬ 
tial  TPAL.  Using  Figure  19  we  arrive  at  the  laser  rate  equations  listed  in  Equation  30. 
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Quenching  rates  are  not  a  significant  process  since  the  pump  level  (52D5/2)  is  also  the 
upper  laser  level  to  the  5.2  firn  laser.  This  laser  then  populates  the  upper  laser  level 
of  the  blue  laser. 


Figure  19.  The  3-level  TPAL  system.  The  pump  laser  need  only  pump  to  the  virtual 
state,  as  second  order  absorption  sends  population  to  the  52D5/2  state.  An  IR  laser 
then  sends  population  to  the  upper  laser  level  of  the  blue  TPAL. 


dns 

dt 

dn2 

dt 


<j> 


ni  -  T IRn3 


^IRn  3  —  ^21^2 


(30) 

(31) 


where  is  the  branching  ratio  for  the  first  laser  transition  and  A21  is  the  spon¬ 
taneous  emission  for  the  upper  laser  level.  Applying  a  steady  state  condition  to 
Equation  30,  we  can  define  an  approximation  to  the  laser  threshold  condition  to  be, 


Assuming  that  we  have  no  losses  in  the  system,  a  — >  0,  we  can  define  the 
population  ratio  in  Equation  32  to  be  2.  Equation  32  now  shows  the  importance 
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of  knowing  the  TPA  cross-section  as  a  first  step  in  developing  a  TPAL.  For  a  TPA 
cross-section  of  6.75  x  1CT21  the  calculated  threshold  intensity  is  Ith  =  26.0 
This  is  about  26  times  higher  than  a  typical  DPAL  value  for  threshold.  This  value 
indicates  a  difficulty  in  developing  a  TPAL  laser  at  higher  alkali  concentrations  for  a 
cw  pump  source.  However,  a  pulsed  dye  laser  has  sufficient  pump  intensities  to  realize 
a  TPAL  laser.  Under  typical  laser  operations  Doppler  widths  must  be  accounted  for, 
which  is  why  our  value  for  the  TPA  cross-section  is  the  most  practical. 
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IV.  XPAL 


4.1  Introduction 

The  Diode  Pumped  Alkali  Laser  (DPAL)  and  its  cousin,  the  Exciplex  Pumped 
Alkali  Laser  (XPAL),  have  received  considerable  attention  for  high  power  applications. 
The  DPAL  system  is  optically  pumped  in  the  core  of  the  D2  absorption  line  of  K, 
Rb,  or  Cs  and  lases  on  the  Di  line  after  collisional  fine  structure  mixing  between  the 
excited  2p3/2,i/2  states  [17,  28,  56,  65].  In  contrast,  the  XPAL  system  is  pumped  in 
the  blue  satellite  wing  of  the  D2  line  and  lases  either  on  the  D2  or  Di  line  [34,  35,  36, 
37,  38].  A  pulsed  XPAL  was  first  demonstrated  in  2009  by  pumping  atomic  cesium 
in  a  mixture  of  500  Torr  of  argon  to  induce  the  satellite  blue  peak  near  837  nm 
and  100  Torr  of  ethane  to  collisionally  relax  the  dissociated  Cs  2P^,/2  population  to 
the  upper  laser  level,  2P\/2  [35].  Lasing  on  the  2P\/2  — >2  S i/2  D\  line  at  894.3  nm 
is  achieved  across  a  broad  band,  ~  3  nm  of  pump  wavelengths,  in  this  five-level 
version  on  the  XPAL  system.  Threshold  was  observed  at  40  —  120  J/ pulse  and  slope 
efficiencies  of  up  to  5%  were  achieved,  depending  on  cesium  concentration  [34,  35,  38]. 
More  recently,  a  four-level  version  of  the  XPAL  was  demonstrated  with  high  quantum 
efficiency,  >  98%,  and  no  requirement  for  hydrocarbons  [36,  37].  The  four- level  Cs-Ar 
system  exhibits  a  threshold  at  130—250  J  in  a  4  ns  pulse  with  an  efficiency  of  about  1% 
and  has  been  scaled  in  pump  rate  to  2  — 3  times  threshold  [4,  36,  37].  Lasers  operating 
at  852.1  nm  in  alkali-rare  gas  mixtures  of  Cs-Ar,  Cs-Ar/Xe,  and  Rb-Ar  at  pressures  of 
500  Torr  have  also  been  demonstrated  at  higher  pump  energy,  <  40  mJ,  and  longer 
pulse  durations,  6  —  11  ns  [33].  The  threshold  is  higher,  1.5  —  12  mJ  /pulse,  and  slope 
efficiency  less,  ~  0.4%,  presumably  due  to  a  larger  pump  spot  size.  In  the  present 
work,  we  seek  to  extend  the  four-level  Rb-Ar  XPAL  performance  characterization  to 
higher  pressures  (2500  Torr),  scale  the  pump  rate  to  >  8  times  threshold,  investigate 
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the  laser  dynamics,  compare  static  gain  cells  with  heat  pipes,  and  demonstrate  a 
Rb-Kr  XPAL  variant. 

The  potential  energy  surfaces  for  several  alkali  metal  vapor-rare  gas  pairs  have 
recently  been  computed  at  the  spin-orbit  multi-reference  configuration  interaction 
level  [46].  The  Rb-Ar  surfaces  for  the  lowest  electronic  states  are  illustrated  in  Fig¬ 
ure  20.  In  the  separated  atom  limit  the  X2EA2  surface  corresponds  to  the  2Si/2 
ground  state  energy  level  of  the  alkali,  the  A2 IR/2  surface  corresponds  to  the  2P\/2 
energy  level,  and  the  A2II3//2  and  B2 E]y2  surfaces  correspond  to  the  2P3/2  energy  level. 
The  X2E+2  surface  is  mostly  repulsive,  with  a  shallow  well  (49  crn~l )  in  the  region 
of  R  fa  0.58  nm.  The  B2 E]y2  surface  also  exhibits  a  shallow  well  (7.3  crrT1)  in  the 
asymptotic  limit  of  R  ~  10.5  nm,  as  well  as  a  shoulder  on  the  repulsive  wall  at  smaller 
values  of  R.  Dipole  allowed  transitions  of  the  Rb-Ar  collision  pairs  from  the  ground 
state  to  the  dissociative  B2 E]y2  surface  produces  a  satellite  on  the  blue  side  of  the  D2 
atomic  line.  The  blue  wing  for  the  Rb-Ar  system  has  been  experimentally  observed 
at  755.5  or  754  nm  [18].  Lasing  emission  on  the  Rb  D2  line  occurs  at  780.2  nm. 


Internuclear  Distance,  R  (A) 

Figure  20.  The  first  four  Rb-Ar  surfaces  and  three  lowest  atomic  Rb  states  [46]. 
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Most  DPAL  and  all  XPAL  lasers  have  previously  employed  static  gas  cell  at  a  fixed, 
near  atmospheric  pressure  for  the  gain  medium.  Heat  pipe  configurations  offer  higher 
alkali  concentrations,  variable  buffer  gas  pressures,  and  avoid  alkali  condensation  on 
cell  windows  [30].  A  3  atm,  hydrocarbon  free  heat  pipe  was  used  to  scale  pulsed  pump 
intensities  in  a  K-He  DPAL  to  greater  than  10  MW/ cm2  while  retaining  greater  than 
10%  slope  efficiency  [31].  Stimulated  Raman  and  multi-photon  excitation  processes 
have  been  observed  at  these  very  high  pump  intensity  conditions  [42],  The  utility  of 
the  heat  pipe  configurations  for  XPAL  devices  are  explored  in  the  present  work. 

4.2  Experimental 

A  schematic  diagram  of  the  experimental  set-up  is  shown  in  Figure  21a.  LInlike 
previous  studies  demonstrating  the  XPAL  system  [4,  33,  34,  35,  36,  37,  38],  this  study 
used  a  glass  cell  connected  to  a  gas  handling  system  or  a  heat  pipe  [30]  to  evaluate 
high  pressure  conditions.  In  previous  studies,  a  sealed  cell  with  a  fixed  amount  of 
rare  gas  and  a  few  micro-grams  of  alkali  was  employed  [4,  33,  34,  35,  36,  37,  38].  The 
glass  cell  used  in  the  current  experiment  has  a  vertical  “finger”  located  on  the  bottom 
of  the  cell  and  centered  between  the  windows  with  about  a  gram  of  Rb.  Heat  tape 
was  wrapped  around  this  finger,  monitored  with  a  resistance  temperature  detector 
and  maintained  to  1  °C  by  an  Omega  temperature  controller.  The  vapor  pressure  of 
Rb  at  T  =  453  —  513  K  yields  concentrations  of  0.34  —  3.3  x  1015  ato™s ,  The  17  cm 
long  glass  cell  rests  in  a  stainless  steel  oven  at  5  — 10  °C  above  the  finger  temperature 
to  prevent  window  condensation.  Vacuum  («  1  Torr )  and  buffer  gas  handling  was 
controlled  through  a  valve  and  glass-to-metal  interface  on  the  top  the  cell. 

The  stainless  steel  heat  pipe  illustrated  in  Figure  21c  was  substituted  for  the  glass 
cell  for  several  experiments  [30].  About  2  g  of  Rb  was  loaded  using  a  nitrogen  purged 
glove  box.  A  stainless  steel  oven  is  clamped  around  the  center  of  the  heat  pipe.  Two 
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Figure  21.  The  Rb-Ar  apparatus.  Panel  (a)  display  the  optical  path  and  important 
elements  of  the  beam  train.  Panel  (b)  describes  the  geometry  of  the  glass  test  cell  used 
for  part  of  the  experiments.  Panel  (c)  explains  the  geometry  and  important  features 
of  the  heat  pip  configuration  used  for  the  majority  of  the  results. 


water-cooled  (20  °C)  jackets  towards  the  ends  of  the  pipe  condense  the  alkali  and  a 
wire  mesh  inside  the  heat  pipe  wicks  the  liquid  alkali  back  to  the  center.  Although 
the  heat  pipe  is  38  cm  long,  the  active  gain  medium  between  the  cooling  jackets 
measures  about  24  cm  in  length.  Heat  pipe  temperature  of  398  —  448  K  yields  alkali 
concentrations  of  0.22  —  2.76  x  1014  atcc^ls ■  The  rare  gas  purity  was  99.998%  and 
the  maximum  pressure  was  2500  Torr.  Pressure  was  measured  using  an  MKS  model 
626A13TAE  capacitance  manometer.  The  alkali  used  was  natural  isotopic  abundance 
rubidium. 

A  Q-switched  Quanta  Ray  Nd:YAG  laser  pumped  a  Sirali  model  PRSCD18  pulsed 
dye  laser  to  produce  up  to  23.5  mJ  in  a  8.3  — 9.7  ns  pulse  at  10  Hz.  Exciton  LDS  765 
dye  in  propylene  carbonate  and  methanol  provided  a  tuning  range  of  735  —  795  nm, 
with  a  peak  at  about  762  nm.  The  spot  size  was  determined  by  capturing  burn 
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images  of  the  incident  pump  pulse  prior  to  entering  the  test  cell.  For  energies  ranging 
from  0.35  —  21.2  mJ/pulse  spot  sizes  ranged  from  0.13  —  0.20  cm 2 .  A  fast  (1  ns) 
Si  photodetector  and  3  GHz  oscilloscope  monitored  pump  and  XPAL  laser  pulses. 
A  780  nm  notch  filter,  with  a  10  nm  bandwidth,  was  used  to  filter  out  the  pump 
wavelength.  An  Ocean  Optics  spectrometer  was  used  to  analyze  the  pump  and  XPAL 
wavelengths.  Average  output  power  was  recorded  on  a  Coherent  FieldMax-TO  power 
meter. 

The  pump  beam  enters  the  XPAL  cavity  by  way  of  a  Polarizing  Beam  Splitter. 
The  cavity  mirrors  consisted  of  a  25  cm  high  reflector  and  a  flat,  30%  transmission 
output  coupler,  spaced  by  53  cm.  The  cavity  round  trip  time  is  3.5  ns.  Sapphire 
windows  at  normal  incidence  were  used  on  the  ends  of  the  heat  pipe.  Insertion  of 
the  heat  pipe  and  beam  splitter  into  the  cavity  introduces  a  large  loss  of  «  60  —  70% 
per  pass  for  a  HeNe  632.8  nm  probe.  The  heat  pipe  loss  increases  significantly  at 
pressures  above  1500  Torr,  particularly  at  high  temperature  (448  K),  as  discussed 
below. 

4.3  Results 

The  current  laser  demonstrations  were  conducted  using  the  four  level  variant  of 
the  XPAL  system.  Figure  22  illustrates  the  spectrum  of  the  XPAL  power  using  the 
heat  pipe  at  500  Torr  of  argon  at  a  cell  temperature  of  494  K .  The  pump  source  is 
tuned  to  755.5  nm  in  the  far  blue  wing  of  the  D2  line  and  significant  pump  energy  is 
transmitted  through  the  cell.  Strong  lasing  is  observed  on  the  Rb  52P3/2  — *  52<Si/2 
transition  near  780  nm.  No  laser  emission  is  observed  on  the  52Pi/2  — »  52S'i/2 
transition  near  795  nm.  A  780  nm  notch  filter,  with  a  10  nm  bandwidth,  was  used 
to  filter  out  the  pump  wavelength  in  subsequent  measurements. 

The  dependence  of  the  XPAL  output  energy  on  the  pump  wavelength  is  shown 
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Figure  22.  An  Ocean  Optics  Spectrometer  was  used  to  verify  emission  of  the  laser 
level. 

in  the  excitation  spectra  of  Figure  23.  Despite  the  ~  20%  pulse-to-pulse  variation 
in  both  pump  and  XPAL  output  energy,  the  peaks  at  755.0  ±  0.5  nm  for  the  Rb-Ar 
system  and  759±0.3  nm  for  Rb-Kr  are  clearly  evident.  Lasing  threshold  for  the  Rb-Ar 
system  is  observed  on  the  blue  side  at  752.3  nm  where  the  absorption  cross-section  is 
reduced  by  a  factor  of  3.3  —  4.6  from  the  peak,  depending  on  temperature  [18].  Lasing 
is  sustained  for  pump  wavelengths  longer  than  760  nm  where  the  absorption  reaches 
a  minimum  before  extending  into  the  core  of  the  D2  line.  In  contrast,  the  lower 
pressure  Rb-Kr  system  supports  lasing  in  a  narrow  range  of  757.1  —  760.4  nm.  The 
reduced  performance  for  the  Rb-Kr  system  in  Figure  23  reflects  the  lower  pressure 
conditions  rather  than  the  relative  merit  of  the  collision  partner.  Comparing  Ar  and 
Kr  performance  at  common  500  Torr  and  398  K  conditions,  yields  output  energies 
of  0.92  and  3.2  J,  respectively.  This  is  qualitatively  consistent  with  the  observed 
absorption  cross-section  for  Kr  being  about  70%  greater  than  for  Rb-Ar  [18]. 

The  scaling  of  output  energy  with  absorbed  pump  energy  for  the  heat  pipe  con¬ 
figuration  is  explored  using  two  approaches  in  Figure  24.  First,  the  energy  per  pulse 
as  measured  by  the  power  meter  and  inferred  from  the  10  Hz  repetition  rate  is 
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Figure  23.  Excitation  spectra  of  both  the  Rb-Ar  and  Rb-Kr  pair.  These  peaks  centered 
at  755.0 ±0.5  nm  and  759  ±0.3  nm  respectively.  The  Rb-Ar  laser  excitation  curve  is  taller 
and  wider  than  the  Rb-Kr  curve  because  of  the  conditions  (temperature  and  pressure) 
by  which  they  were  observed. 


measured  for  both  the  pump  and  XPAL  lasers.  The  fraction  of  the  incident  pump 
energy  absorbed  by  the  Rb  vapor  is  determined  by  comparing  the  transmitted  power 
at  the  peak  of  the  blue  wing  feature  at  A  =  755  nm  and  the  off-resonance  at 
A  =  750  nm.  For  1500  Torr  of  argon  and  T  =  423  K  the  fraction  absorbed 
ranges  from  0.804  —  3.24  mJ  /pulse  as  the  incident  pump  energy  is  scaled  from 
5.4  —  21.8  mJ /pulse.  A  linear  fit  to  this  data  is  shown  in  Figure  24  and  yields  a 
threshold  of  0.38  ±  0.02  mJ /pulse  of  absorbed  energy  or  2.55  mJ /pulse  incident  en¬ 
ergy.  For  the  10  ns  pulse  focused  to  an  area  of  0.01  cm2,  this  corresponds  to  an 
average  pump  intensity  of  25.5  ^r.  The  fit  slope  efficiency  is  0.236  ±  0.003  %. 

A  second  approach  to  developing  the  results  of  Figure  24  employs  the  full  tem¬ 
poral  dependence  of  the  pump  and  XPAL  output  pulses.  Using  the  fast  photodiode 
and  normalizing  the  temporally  integrated  pulse  area  to  the  observed  average  energy 
provides  the  instantaneous  intensity  curves  in  Figure  24.  For  an  incident  pump  en- 
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Figure  24.  Two  methods  for  obtaining  the  slope  efficiency  are  shown.  The  data  with 
the  linear  fit  was  collected  by  the  common  method  of  measuring  the  absorbed  pump 
energy  and  output  laser  energy  at  certain  intervals  of  pump  intensity.  The  more  non- 
traditional  method  is  to  use  a  single  laser  pulse  to  sketch  out  all  needed  energies.  A 
photon  build-up  time  was  determined  from  the  geometry  of  the  set-up.  When  factored 
out,  the  hysteresis  collapses  to  a  steady  state  condition.  The  discrepancy  in  intercepts 
is  probably  due  to  the  different  methods  in  collecting  the  data.  Only  the  slope  efficiency 
was  collected  from  this  method  in  this  study. 

ergy  of  21.8  mJ/ pulse ,  the  instantaneous  intensity  traces  an  oval  shape  propagating 
counter-clockwise  as  time  progresses.  The  full  temporal  dynamics  of  the  two  pulses 
can  be  organized  by  a  simple  temporal  delay  between  the  pump  and  XPAL  pulses  of 
about  Tb  =  3  ns,  which  can  be  attributed  to  the  photon  buildup  time.  Introducing 
this  delay  reduces  the  hysteresis  to  a  single  linear  relationship,  consistent  with  the 
average  energy  measurements.  A  pump  pulse  with  lower  energy  of  12.8  mJ /pulse 
yields  similar  results.  Thus,  the  full  scaling  relationship  can  be  interrogated  from  a 
single,  high  energy  pulse.  It  is  worth  noting  that  some  nonlinear  behavior  is  observed 
near  threshold  in  all  the  plots  of  Figure  24,  consistent  with  the  upward  curvature 
(increasing  slope  efficiency)  observed  in  the  earlier  XPAL  demonstrations  [34,  35,  38]. 

The  photo  build-up  times  observed  in  the  glass  cell  for  all  the  argon  pressures  and 
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cell  temperatures  are  provided  in  Figure  25.  The  delay  decreases  from  3.7  to  2.8  ns 
as  the  Rb  concentration  increases  by  a  factor  of  9.7  (achieved  by  changing  the  cell 
temperature  from  453  —  513  K)  The  photon  build-up  time  is  inversely  proportional  to 
the  gain,  and  the  XPAL  inversion  depends  on  both  the  rubidium  concentration  and 
the  argon  pressure  (due  to  the  increased  blue  wing  absorption)  [4],  Indeed,  Figure  25 
demonstrates  an  approximately  linear  relationship.  The  gain  is  clearly  large  and  is 
estimated  as  7  ~  =  18  —  24—,  where  the  factor  of  20  is  a  typical  value  for  well 

saturated  lasers.  Delays  using  the  heat  pipe  configuration  at  higher  argon  pressures 
(less  than  2500  Torr )  and  longer  active  region  are  even  less,  1.0  —  1.9  ns. 


Figure  25.  The  photon  build-up  time  in  the  glass  cell  configuration  was  shown  to 
decrease  with  increasing  alkali  concentration  as  expected. 

The  observed  slope  efficiencies  increase  as  the  Rb  and  Ar  concentrations  increase, 
as  shown  in  Figure  26.  The  absorption  of  the  pump  photons  in  the  blue  wing  of  the 
D2  line  depends  on  both  the  alkali  and  rare  gas  concentrations  [4,  35,  37,  38,  48]. 
An  analytic  model  for  the  standard  DPAL  system,  predicts  a  slope  efficiency  which 
depends  exponentially  on  the  absorbance,  A:  [28] 

Vs  =  7/o(l  -  e'A)  (33) 
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A  fit  of  Equation  33  to  the  data  in  Figure  26  yeilds  an  effective  absorption  cross- 
section  of  A  =  a  n  L  of  o  =  8.16  ±  1.8  x  ICR17  cm 2  at  argon  pressures  of  500  and 
700  Torr ,  which  is  approximately  linearly  dependent  on  buffer  gas  pressure.  This  is 
about  10-4  of  the  peak  D2  cross-section,  consistent  with  prior  asymmetric  lineshape 
observations  [33,  35,  37,  38,  48].  The  efficiency  in  the  heat  pipe  configuration  is 
dramatically  higher,  despite  the  lower  Rb  concentration.  However,  at  high  pressure 
and  temperature  efficiency  becomes  quite  poor,  as  discussed  below.  These  lower 
efficiency  conditions  are  not  included  in  the  plot  of  Figure  26.  The  current  Rb-Ar 
slope  efficiencies  are  lower  than  previously  reported  for  the  Cs-Ar  system  of  fa  1%  [4], 
consistent  with  the  qualitative  observations  reported  by  Reaclle  [35,  36,  37,  38]. 


Figure  26.  We  extract  an  effective  absorption  cross-section  by  fitting  the  function 
derived  by  Hager.  [28] 


The  number  of  pump  photons  per  pulse  is  high,  8.29x  1016  photons  at  21.8  mJ /pulse. 
The  number  of  available  Rb  atoms  in  the  pumped  volume  of  the  glass  cell  is  lower, 
2.22  x  1013  to  2.76  x  1014  atoms,  depending  on  the  cell  temperature.  The  number 
of  Rb-Ar  pairs  is  lower  still  and  can  be  estimated  from  the  potential  surfaces  [46]. 
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This  effect  is  accounted  for  when  considering  the  reduced  cross-section  when  pump¬ 
ing  in  the  blue  wing.  The  fraction  of  the  incident  photons  absorbed  is  low  and  the 
optical-to-optical  efficiency  suffers  from  the  highly  bleached  conditions. 

We  anticipated  improved  performance  as  both  Rb  and  Ar  density  increased  further 
and  the  heat  pipe  provided  capability  to  extend  the  range  to  T  <  800  K  and  P  < 
2500  Torr.  The  slope  efficiencies  observed  for  the  heat  pipe  configuration  are  provided 
in  Figure  27.  Surprisingly,  no  lasing  was  achieved  in  this  configuration  for  T  >  448  K . 
Slope  efficiencies  increase  initially  with  pressure  and  temperature,  but  degrade  rapidly 
above  428  K  and  2000  Torr. 


Figure  27.  Slope  efficiencies  for  the  various  test  conditions  performed  with  the  heat 
pipe.  An  obvious  degradation  in  performance  is  shown  for  high  pressures  and  temper¬ 
atures.  We  don’t  believe  this  is  a  real  XPAL  effect,  but  rather  a  consequence  of  not 
fully  understanding  the  heat  pipe  operational  conditions. 


We  visually  observed  condensation  of  alkali  vapor  and  the  production  of  suspended 
aerosols  in  the  heat  pipe  at  the  higher  pressures  and  temperatures,  only  when  the 
pump  laser  was  incident  on  the  medium.  Figure  28  illustrates  the  optical  transmission 
for  a  low  power  HeNe  laser  (632.8  nm)  as  a  function  of  HeNe  pump  energy  for  several 
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heat  pipe  conditions.  The  transmission  is  high  and  independent  of  HeNe  pump  energy 
for  the  lower  pressures  and  temperatures.  However,  the  transmission  is  very  poor 
at  elevated  pressure  and  temperature  for  incident  pump  energies  exceeding  0.5  nJ . 
Indeed,  there  is  an  abrupt,  unstable  transition  to  low  transmission  when  the  heat 
pipe  is  illuminated  for  more  than  5  seconds,  as  illustrated  in  Figure  29.  Apparently, 
the  heat  load  induced  by  the  pump  beam  induces  a  phase  change  in  the  heat  pipe  at 
higher  operating  pressures  and  temperatures,  limiting  further  scaling. 


Pump  Energy,  Ep  (pJ) 


Figure  28.  Transmission  of  a  low  power  HeNe  laser  down  the  optical  axis  of  the  heat 
pipe  as  a  function  of  pump  energy,  alkali  concentration  and  pressure.  Each  curve  is 
labeled  with  the  condition  the  data  was  collected  under. 


Finally,  a  Rb-Kr  XPAL  system  was  demonstrated  using  the  heat  pipe  configu¬ 
ration.  The  instantaneous  slope  efficiencies  for  the  Kr  buffer  gas  are  provided  in 
Figure  30.  Unfortunately,  no  lasing  was  observed  for  the  higher  pressure  and  temper¬ 
ature  conditions.  However,  the  production  of  aerosols  for  Kr  system  appears  to  occur 
more  readily  and  is  likely  the  explanation  for  the  degraded  performance.  Indeed, 
the  Rb-Kr  system  exhibits  improved  slope  efficiencies  over  the  Rb-Ar  system  for  the 
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Figure  29.  Transmission  of  a  low  power  HeNe  laser  down  the  optical  axis  of  the  heat 
pipe  as  a  function  of  time.  This  shows  the  degrading  mechanism  to  have  a  time  depen¬ 
dent  function  that  increases  with  pressure. 


lowest  pressures  and  temperatures. 


4.4  Conclusion 

The  performance  characterization  of  the  exciplex  pumped  alkali  laser,  XPAL, 
has  been  extended  to  pressure  above  1  atmosphere  and  pump  intensities  exceeding 
15  in  the  present  work.  Threshold  is  observed  at  about  0.5  mJ  absorbed  energy 
(3  mJ /pulse  incident  energy)  and  the  system  has  been  scaled  to  greater  than  8x 
threshold.  The  laser  dynamics  are  consistent  with  quasi-cw  behavior,  requiring  only 
a  photon  build-up  time  of  1  —  3  ns  to  explain  the  temporal  nature  of  the  output. 
However,  a  cw  XPAL  demonstration  remains  elusive,  due  to  the  high  pump  intensities 
required.  Extending  performance  to  higher  pressures  and  temperatures  using  a  heat 
pipe  will  require  control  of  condensation  induced  by  laser  heating.  The  first  XPAL 
system  using  Rb-Kr  pairs  was  presented  with  performance  potentially  exceeding  the 
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Figure  30.  Slope  efficiencies  for  the  Rb-Kr  XPAL. 


Rb-Ar  system. 


58 


V.  Conclusion 


This  work  has  demonstrated  for  the  first  time  a  measurement  of  the  Doppler- 
broadened  TPA  cross-section  for  the  52S'i/2  — >  5 2D5/2  transition  in  naturally  occur¬ 
ring  Rb  using  direct  absorption  methods.  Our  reported  statistical  average  hyperhne 
broadened  TPA  cross-section  is  7.6  x  ICR21  This  work  also  provides  an  intensity 
model  that  accounts  for  the  Gaussian  propagation  term  and  the  TPA  cross-section 
absorbing  term,  which  allows  for  absorbance  values  to  be  converted  into  TPA  cross- 
section  units  of  Although  the  values  calculated  by  Zapka  et  al.  and  Saha  et 
al.  were  5  and  6.5  times  our  calculated  value  respectively,  we  believe  that  this  dis¬ 
crepancy  is  largely  accounted  for  in  the  lineshape.  Both  of  their  experiments  used  a 
sub-Doppler  technique,  whereas  our  direct  absorption  method  and  Voigt  model  used 
the  Doppler  widths  at  high  temperatures  to  generate  sufficient  alkali  concentrations. 
The  values  calculated  by  Marinescu  et  ah  and  Collins,  however,  were  orders  of  mag¬ 
nitude  different  and  it  is  believed  they  incorporated  subtle  errors  into  their  work.  The 
calculated  TPA  cross-section  in  this  work  indicates  a  required  pump  intensity  thresh¬ 
old  that  is  approximately  25  times  greater  than  a  typical  DPAL  system.  Although 
these  intensities  are  achievable  using  pulsed  lasers,  the  likelihood  of  developing  a  cw 
TPAL  system  is  greatly  decreased.  The  value  calculated  in  this  work  becomes  the 
most  practical  reported  number  when  considering  the  feasibility  of  a  TPAL  system 
because  the  lineshape  dependence  of  the  cross-section.  In  typical  laser  operations  the 
Doppler  width  must  be  accounted  for. 

In  the  area  of  XPALs,  this  work  was  the  first  to  develop  a  working  Rb-Ar  and 
Rb-Kr  XPAL  in  a  heat  pipe  configuration.  These  particular  XPALs  optically  pumped 
the  blue  satellite  of  the  D2  line  from  the  ground  Rb-Ar  or  Rb-Kr  collision  pair  to  the 
dissociative,  R2S^2  state  produce  laser  emission  at  780.2  nm.  We  achieved  lasing 
for  a  band  of  pump  wavelengths  of  752.6  —  760.0  nm  for  the  Rb-Ar  system  and 
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757.7  —  759.8  nvn  for  the  Rb-Kr  system.  This  work  was  also  the  first  to  characterize 
XPAL  laser  operation  at  pressures  above  1  atmosphere  and  pump  intensities  exceeding 
15  Threshold  was  observed  at  about  0.5  mJ  absorbed  energy  (3  mJ/pulse 

incident  energy)  and  the  system  has  been  scaled  to  greater  than  8  times  threshold. 
This  is  significantly  less  than  the  values  calculated  for  the  Cs-Ar  systems  by  Readle 
et  ah,  which  ranged  from  40  pJ  at  444  K  to  140  pJ  at  468  K.  Our  operating 
temperatures  and  pressures  were  limited  by  a  condensation  of  the  alkali  due  to  the 
heat  pipe  configuration.  This  effect  was  exacerbated  by  laser  induced  heating  and 
increased  alkali  vapor  pressures.  The  laser  dynamics  are  consistent  with  quasi-cw 
behavior,  requiring  only  a  photon  build-up  time  of  1  —  3  ns  to  explain  the  temporal 
nature  of  the  output.  However,  a  cw  XPAL  demonstration  remains  elusive,  due  to 
the  high  pump  intensities  required.  Extending  performance  to  higher  pressures  and 
temperatures  using  a  heat  pipe  will  require  control  of  condensation  induced  by  laser 
heating.  Finally,  slope  efficiencies  increased  with  both  Rb  and  Ar  concentrations  and 
exceed  0.25%  using  the  heat  pipe  configuration. 

5.1  Contributions 

The  TPA  cross-section  leads  to  efforts  in  developing  an  efficient  TPAL  system. 
This  laser  system  has  the  potential  to  benefit  Naval  communications  for  submarines 
in  shallow  water.  The  420  nm  wavelength  is  able  to  penetrate  further  than  most 
optical  wavelengths  through  the  water  absorbing  bands.  Pumping  to  other  two- 
photon  transitions  in  Rb  will  allow  access  to  other  micron-range  laser  systems,  which 
could  also  benefit  advanced  IR  countermeasure  technology. 

The  XPAL  research  provides  important  data  to  the  growing  laser  community 
seeking  for  the  next  non-chemical  high-power  laser  system.  The  results  from  this 
work  suggests  exciting  areas  to  explore  using  the  heat  pipe  configuration,  and  provides 


60 


lessons  learned  from  the  unexpected  results  presented. 


5.2  Future  Work 

Many  areas  have  been  identified  as  potential  research  areas  throughout  the  course 
of  this  research  in  the  TPA  cross-section.  First  of  all,  it  was  the  intention  of  this  re¬ 
search  to  calculate  the  pressure  broadening  rates  by  extracting  Lorenztian  linewidths 
from  a  series  of  pressure  broadened  TPA  transitions  in  Rb.  This  will  greatly  enhance 
the  understanding  of  the  contributing  factors  of  a  lineshape  resulting  from  a  TPA 
process.  The  next  obvious  step,  after  calculating  the  broadening  rates,  is  to  build 
a  TPAL  cavity  and  test  various  designs  while  monitoring  and  modeling  the  perfor¬ 
mance.  The  most  interesting  aspect  of  this  transition  is  the  simultaneous  emission  of 
two  laser  wavelengths.  This  fact  alone  could  drive  varying  ideas  of  cavity  configura¬ 
tion.  Glass  cells  were  used  in  the  measuring  of  the  TPA  cross-section,  but  considering 
the  use  of  a  heat-pipe  is  also  very  promising.  A  heat-pipe  configuration  would  not 
experience  the  degradation  in  performance  at  high  temperatures  that  this  research 
encountered.  Using  a  heat-pipe,  one  could  even  redo  much  of  this  work  and  then 
extend  to  higher  temperatures  and  pressures. 

For  the  XPAL  work,  demonstrating  and  reporting  the  adequate  operating  condi¬ 
tions  for  the  heat  pipe  configuration  would  be  very  useful  for  analyzing  high  pressure 
XPAL  performance.  Avoiding  large  temperature  gradients  is  a  challenge  when  work¬ 
ing  at  high  temperatures  in  a  heat  pipe.  Understanding  the  proper  setting  of  the 
cooling  jackets  may  have  a  significant  contribution.  Once  the  heat-pipe  operation  is 
better  understood,  it  is  recommended  that  various  alkali-rare  gas  combinations  are 
tested  and  analyzed  for  performance. 
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Appendix  A.  Developing  the  General  Equation  and  the 

Numerical  Solution 

The  intensity  modeled  by  Equation  19  is  the  result  of  developing  a  relationship 
that  accounts  for  the  two  competing  processes  found  within  the  alkali  test  cell.  The 
first  process  describes  the  intensity  increasing  due  to  the  decreasing  beam  waist, 
and  the  second  describes  the  intensity  attenuated  through  the  absorption  medium. 
Figure  31  illustrates  the  processes  needed  to  develop  a  model  for  the  intensity  profile 
in  the  test  cell.  As  the  shaded  region  travels  in  the  z-direction,  the  intensity  increases 
linearly  in  /  as  it  approaches  the  beam  waist  and  decreases  quadratically  in  /  as  the 
flux  of  photons  increases  and  invokes  stronger  TPA. 

A 


L 


Figure  31.  The  laser  beam  passes  through  a  focus  as  it  travels  through  the  Rb  Vapor. 
An  accurate  calculation  of  the  TPA  cross-section  requires  a  model  of  I(z). 

The  first  process  is  a  result  of  a  Gaussian  beam  propagating  through  a  focus  and 
the  second  arises  from  the  two-photon  version  of  the  Beer-Lambert  Law.  The  general 
equation  describing  the  behavior  of  these  processes  can  be  written  as, 

I(z  +  dz )  =  I(z)  —  a2(j>NI2{z)dz  +  F(z)I(z)dz.  (34) 
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I (z)  is  a  known  intensity  at  some  point  z  and  I(z  +  dz )  is  the  intensity  at  a  later 
point  dz  away.  The  second  term  is  the  TPA  model  and  the  third  term  is  the  focusing 
behavior  of  the  Gaussian  beam  propagating  through  free  space,  which  was  derived 
from  Equation  11.  It  is  necessary  to  subtract  a  factor  of  h  from  the  propagation 
distance,  z,  in  the  Gaussian  term  in  order  to  linearly  add  the  two-photon  term  to  the 
Gaussian  term.  Using  Equations  11  and  34,  along  with  careful  algebra,  we  arrive  at 
the  following  equation  which  this  work  refers  to  as  the  “Governing  Equation” . 

d~ir  =  + 

where, 

cu(L)=co20  [1  +  a2{L-H)2] 

and, 

doj(L)  coo  a2(L  —  H) 

dL  ~  [l  +  a2{L- H)2}1/2 

Substituting  the  relations  of  Equations  36  and  37  into  Equation  35,  yields  the 

differential  equation  that  is  a  function  of  a  and  L. 


(36) 


(37) 


2  du(L  ')  f  r 


u(L)  dL  \UJ2(L) 


-  1 


m 


(35) 


dI(L) 

dL 


g(L ) 


a  2 


a24>NI2{L)  + 


2  a2(L-H) 

[1  +  a2(L-H)2  V^o  t1  +  «2(L  -  H)2} 


I{L)  (38) 


U!  =  V(t)  +  s(i)/(i) 


(39) 


We  take  a  moment  to  divert  from  the  full  development  in  order  to  find  a  definition 
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for  I\{L)  from  Equation  17.  We  do  this  by  forcing  a  one-photon  approximation  to 


ai 


Equation  39,  by  setting  a2/2(L)  — >  a2/0  /(T),  such  that 


dI{L) 

dL 


=  -aJ{L)  +  g(L)I(L)  =  (-a,  +  g(L))I(L ) 


(40) 


Now  we  have, 


and, 


dI{L) 


-a\  +  g(L))dL 


dI(L ) 

w 


h(L) 


—  —diL  +  /  g(L)dL+C 


(41) 


(42) 


This  is  how  we  arrive  at  our  definition  for  I\{L)  in  Equation  17.  This  integral  is 

2 

analytic  and  will  be  useful  to  solve  here.  First,  let  us  further  set  fd  =  ^  and  perform 

U!g 

a  u  —  substitution  such  that  u  =  1  +  a2(L  —  H)2 . 


/3du  f  du  n  f  9  ,  f  du 
-  I  —  =  (3  u  du  -  — 


u 


u 


u 


(43) 


,-l 


=  / 3 — —  —  ln(u )  =  —  —  —  ln{u) 


(44) 


Substituting  u  back  into  Equation  44,  we  get  the  final  answer  for  Ii(L). 


Il(L)  =  ~l  +  a2(L  -  ~  ln  (1  +  °2(L  ~  g>2)  (45) 

It  is  interesting  to  note  that  this  integral  is  analytic  and  is  defined  under  a  one- 
photon  approximation.  We  could  have  arrived  at  this  same  definition  by  simply 
carrying  out  the  Bernoulli  Solution,  but  understanding  the  one-photon  approximation 
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Table  5.  Definitions  to  Bernoulli’s  form. 


Bernoulli  Coefficient 

Governing  Equation 

y 

I(L) 

X 

L 

P(x) 

-  g(L) 

Q(x) 

-  CL'2 

n 

2 

is  more  enlightening.  Equation  35  must  be  solved  using  Bernoulli’s  Solution  since  it 
has  the  same  differential  form  as 


£  +  P(x)y  =  Q(x)yn  (46) 

The  next  logical  step  is  to  rearrange  Equation  39  to  match  the  elements  in  Equa¬ 
tion  46. 


dI^-g(L)I(L)  =  -a2I\L)  (47) 

The  relationships  defined  in  Table  5  is  a  result  of  aligning  Equation  47  with 
Equation  46. 

Next  we  use  Bernoulli’s  solution  in  Equation  48  to  solve  for  I(L). 

yl-ne(l-n)  f  p(x)dx  _  Q  _  \ 


/(L)_1e(_1)^_ff(L)di  =  (-1)  /  -a2e{~l) $ ~9(L)dLdL  +  C 


(49) 


Rearranging  terms  and  substituting  the  definition  for  I\{L)  we  now  can  define  /2(h), 


h(L) 


I(L) 


eh  ^  =  a2  e^dL+C 


(50) 
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Now  our  solution  has  a  form  simple  enough  to  solve  for  the  constant  of  integration. 


m 


e/l(i)  =  a2I2(L)  +  C 


(51) 


The  initial  condition  occurs  at  the  start  of  the  test  cell,  or  when  L=0. 


m 


e7‘<°>  =  02/2(0)  +  C 


(52) 


So, 


C  =  -  02/2(0) 


(53) 


Inserting  Equation  53  into  Equation  51,  we  have  the  start  of  the  complete  solution. 


——e1^  _  a2I2(L)  +  -eh(0)  -  a2I2(0 ) 


(54) 


The  next  several  steps  is  simply  algebra  forming  Equation  54  into  the  form  of  Equa¬ 
tion  8. 


7l_e^.  =  o2(/2(/)-/2(0))  +  7Te'1<°) 


(55) 


ML) _ _PM°)  = 

i(l)  ~  m 


e^V'  -  =  a2(I2(L)  -  J2(0)) 


I0e 


-Mo) 


(56) 


I(L) 


0  ML)e-M0)  -i  =  a2  [I2(L)  -  J2(0)]  J0e-/l(0) 


(57) 


I(L) 


oh(L)-Mo)  =  x  +  a2jQ  [/2(L)  _  /2(o)]  e~M0) 


(58) 
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1 


(59) 


HL)_  [i^D-hm  = 

Jo  '  l  +  a2/0[/2(L)-/2(0)]e-Mo) 

And  finally,  we  get  back  Equation  60. 


m  = 


joe[7i(L)-/i(0)] 


(60) 


1  +  a2<l>NI0  [I2{L)  -  /2(0)]  e-A(0) 

There  are  only  a  few  steps  more  until  we  are  able  to  define  the  eta  term,  rj(a,r). 
Since  I\(L)  has  been  solved  in  Equation  45,  it  is  only  a  matter  of  algebra  to  get  the 
following  relation: 


h(L)  —  h(0)  =  —{3 


1  +  a2H2  1  +  a2(L  —  H)2 


+  In 


1  +  a2H2 

[l  +  a2(L-H)2\ 


(61) 


Substituting  Equation  61  into  Equation  60  yields, 


m  = 


( _ 1  +  _ A  pfi\l  +  a'2H2  1 +a2(L-H)2  ) 

\l+a2(L-H)2J  e 


1  +  a2<pNI0  [ I2{L )  -  J2(0)]  (1  +  a2H2)epy 


(62) 


When  we  invoke  the  symmetry  property  of  the  experiment,  H  —  j,  Equation  62 


reduces  to, 


m  = - - - j  -  v  (63) 

1  +  cr^NI»  lh(L)  -  /2(0)]  (1  +  ^f) 

We  can  now  complete  the  solution  by  inserting  the  definition  for  /2(L)  in  Equation  50 
I2{L)  =  I  eh(L)dL  =  j  e-/3 ( i +^-m2 ) _,n(1' +a2(L~H)^dL  (64) 

and 
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I2(L)  =  /  eIl(L)dL  = 


M 


1  +a'2(L-H)'2 


1  +  a2(L-H)2 
Finally,  we  can  express  the  complete  numerical  method  as: 


dL 


(65) 


m  = 


1  +  C2<pN  IqL  — 


g  ^  \  1+a2 (L  —  H)2  ) 

1  +  a2(L  —  H)2 


dL 


a2L2 


1  + 


77^ U 
1+  4 


(66) 


rj(a,r) 


This  eta  term  can  then  be  conceptually  defined  as  a  multiplicative  factor  to  the 
intensity.  Now  that  we  have  a  better  understanding  of  the  Method  used  to  describe  the 
intensity  profile,  we  must  now  use  that  function  to  transform  the  y-axis  in  Figure  11 
to  the  TPA  unit  of  This  is  done  by  integrating  Equation  66  over  all  r  in  order  to 
properly  model  the  transmitted  power,  Ft-  The  difficulty  is  the  integral  that  can  not 
be  solved  analytically.  The  way  we  integrated  this  equation,  was  by  taking  the  ’’sum 
under  the  curve”.  This  operation  was  performed  in  Mathematica  and  was  written  as, 


/OO 

I(L,r)rdr  (67) 

-OO 

In  essence,  the  Mathematica  code  does  the  same  operation  to  the  intensity  profile 
as  Equation  67.  Specifically  however,  the  code  did  the  following.  First  of  all  you  will 
notice  that  there  is  an  additional  (27rr0.001)  term.  The  number  is  the  dr  that  appears 
in  Equation  67  as  well.  Next,  the  N Integrate  command  numerically  integrates  the 
/2(E)  relation  over  a  span  of  2  from  0  to  L.  Now,  the  order  of  operation  is  such  that 
the  Sum  performs  a  calculation  for  every  value  of  r  from  0  to  0.1  microns  in  steps  of 
0.001  microns.  This  step  is  the  dr  that  was  talked  about  above.  To  get  the  integral 
over  all  r,  this  calculation  is  performed  over  a  range  of  r’s  and  is  added  together 

68 


and  multiplied  by  the  27rr0.001  term.  Everytime  a  calculation  is  done  for  an  r  value, 
the  numerical  integral  is  performed  over  the  range  0  to  L  in  z.  Finally,  this  sum 
is  done  for  a  range  of  predefined  cross-sections,  from  0  to  4  x  1CT20  g^-  in  steps  of 
1  x  10~22  Hundreds  of  PowerE’s  are  calculated. 

The  predicted  Absorbance,  Apred,  can  now  be  determined  using  the  two-photon 
variant  of  Beer’s  Law,  Apre(i  =  Po^°erE  —  1.  Where  P0  is  simply  the  power  measured 
just  before  the  test  cell.  A  range  of  cross-sections  were  chosen  such  that  the  computed 
range  of  absorbance  values  would  include  the  largest  absorbance  value  measured  in 
the  experimental  data.  The  range  of  cross-sections  were  then  plotted  against  the 
predicted  absorbance  values  to  generate  the  curve  in  Figure  32. 


ml  -0  — I 

™0.5  — 

!o.o  - 

a=0.5  - . 
-1.0x10  24  J 


^V///*A'y-,.V,yvVv-Vv.vVv. 


Figure  32.  This  plot  provides  the  information  necessary  to  generate  the  function  that 
will  transform  the  experimental  absorbance  values  into  two-photon  cross-section  units. 


The  resulting  function  from  the  fit  to  this  data  provided  the  conversion  necessary 
for  the  y-axis  in  Figure  11  to  be  transformed  into  TPA  cross-section  units.  The  y-axis 
is  now  in  the  proper  units,  so  that  the  Voigt  fit  that  we  have  developed  will  be  able 
to  extract  the  TPA  cross-section  as  the  amplitude  to  the  data.  This  entire  process 
was  performed  87  times.  Each  transformation  curve  generated  from  the  Mathematica 
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code  was  unique  for  each  of  the  87  cases.  Table  6  lists  the  five  parameters  of  the  fit 
function  for  6  different  test  cases  for  comparison. 


Table  6.  Fit  parameters  tx  and  test  conditions  for  6  of  87  test  cases. 


Case 

to(10-2a) 

<74- 

o 

1 

to 

o 

t2(10”2U) 

t3(io-2U) 

t4(io-20) 

N  (io15ato™s') 

P  (mW) 

f  (cm) 

1 

-2.0866 

9.7517 

0.2061 

-1.335 

-3.254 

4.066 

680 

15.24 

2 

-2.7912 

0.1279 

0.2690 

-94.72 

-6.107 

1.554 

1200 

15.24 

3 

-1.3108 

6.5865 

0.1405 

-1.557 

-1.027 

4.066 

930 

15.24 

4 

-1.4442 

4.3683 

5.3007 

20.22 

-35.72 

3.123 

1200 

20.32 

5 

-0.4894 

2.2202 

2.7547 

-212.5 

-846.3 

6.124 

1200 

20.32 

6 

-2.1502 

6.2500 

7.5389 

43.23 

-67.32 

4.199 

700 

20.32 
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Appendix  B.  Development  of  Other  Methods  to  Model  the 

Intensity 

Five  methods  for  modeling  the  intensity  were  explored.  This  research  first  derived 
the  more  complete  model  which  was  described  in  Appendix  A.  However,  to  better 
understand  the  behavior  of  the  two-photon  absorption  process  with  a  focused  beam, 
and  also  to  provide  a  logical  check  to  the  TPA  cross-section  extracted  from  this 
method,  4  other  lower  fidelity  models  were  derived. 

Model  1:  The  One-Photon  Model. 

The  first  model  is  using  the  one-photon  Beers  Law  equation  and  inserting  the 
relation,  cre//  =  to  determine  the  effect  on  the  transmitted  intensity.  Here  <re// 

is  the  effective  ” one-photon”  cross-section  in  units  of  cm2.  This  method  assumes  that 
the  pump  beam  is  collimated  and  is  only  attenuated  by  the  absorption  medium.  This 
is  shown  in  Equation  68  below. 


Method  1  :  I(L)  =  I0eaa*IoNL  (68) 

dl(z)  =  (Ji^N  I  (z)dz  (69) 

I(z)  =  I0e^NL  (70) 

dl(z )  =  a2<pNI2(z)dz  (71) 


A  =  ^--l 
Pt 


71 


(73) 


1 

3 


hv 


«&♦(»*)  —  alm(Vo)  ^  ^  ^ F"  F'  f B  fi: 


(74) 


F'  F"  iso 


Obviously  the  more  familiar  absorbance  equation,  —In  ( ),  is  derived  from  this 


form. 


Model  2:  Two-Photon,  No  Focus  Model. 

The  second  method  uses  the  two-photon  Beers  Law  variant,  which  originates  from 
the  differential  form  in  Equation  10,  and  assumes  a  collimated  beam.  Method  2,  then, 
is  just  the  solution  to  this  differential  equation. 


Method  2  :  I(L) 


_ I_o _ 

1  +  &2<f>IoNL 


The  next  two  methods  include  a  model  for  the  focused  Gaussian  beam. 


(75) 


Model  3:  Two-Photon,  Focused  Beam,  fd  —  0  Model. 

2 

Recall  from  Appendix  A  that  fd  =  A*.  This  model  simplifies  the  complete  numer- 
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ical  model  from  Equation  66  by  solving  for  the  intensity  only  along  the  optical  axis, 
r  =  0,  or  (d  =  0.  When  this  is  done,  the  exponential  terms  disappear  and  the  integral 
in  the  denominator  becomes  analytical. 


/(L)  = - f - - - , - - 

1  +  (TfyNIo  f  1+a2(L_H)2  dL  (l  + 

From  any  table  of  integrals,  one  can  find  the  following  integral, 


(76) 
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1 


-dx  = 


ax 2  +  bx  +  c~~  \/4ac  -  b 2 

Therefore,  the  integral  in  Equation  76  becomes, 


arctan 


2ax  +  6 
\/ 4ac  —  b2 


(77) 


1  +  o2(L  —  77)2 


dL  = 


x 


^4(a2)(o2772  +  1)  -  (-2 o277)2 

2a2L  +  (-2  a2H) 


x  arctan 


^/4(a2)(a2ff2  +  1)  -  (-2a2//)2 


(78) 


After  simplification,  we  have: 


h(L)13  0  =  —  arctan  [a(L  —  77)] 
a 


(79) 


And, 


^2(0)^  0  =  —  arctan  [—0:77]  (80) 

o 

Invoking  the  symmetry  of  the  experiment  (77  =  L/2)  we  finally  have  a  model  for 
Method  3, 


Method  3  :  7(L)/3=0  = 


1  +  &2<t>NI0 


a 


arctan 


a  L 
2 


(1  +  ^) 


(81) 


7?3(a) 

Model  4:  Two-Photon,  Unfocused  Beam,  Average  Intensity  Model. 

A  fourth  method  was  chosen  to  explore  the  idea  that  knowing  the  average  in¬ 
tensity  through  the  cell  would  provide  the  simplest  accurate  answer.  We  start  with 
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Equation  66  from  Appendix  A.  We  start  by  first  verifying  that  the  following  relation 
is  true: 


PT  =  2  7T  /  I(L,r)r  dr 
Jo 


(82) 


Logically,  this  statement  is  true  if  we  can  prove  the  following  is  also  true, 


P0  =  2  n  /( 0, r)r  dr 

Jo 


(83) 


It  follows  that  for  a  focused  Gaussian  beam,  the  intensity  behaves  as, 


1(0,  r)  =  P°e  ^5) 

7raA(0) 


(84) 


Therefore, 


Pn  =  2  7T 


Pn 


7ro;2(0) 


fOO  r2 

e  “2(°)  r  dr 


(85) 


From  integral  tables  we  find  that, 


■  rfm+lf 

xme~ax  dx=  V  2  ’ 


2a  (^) 


(86) 


Such  that, 


Pn  =  2vr 


Pn 


7to;2(0) 


oj2(0) 


(87) 


And  finally, 


Po  =  Po 


(88) 
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Returning  to  Equation  82,  we  insert  the  relation  we  have  for  I(L,r )  from  Equation  66 
to  get, 


Pt  =  2  7T 


1(0,  r)  r  dr 


0  1  +  a^NL  7(0, r)M^Lj 

Ju 


(89) 


F 


r°°  i 

Pt  =  2tt  /  — —  7(0,  r)  r  dr 

Jo  1  +  7 


(90) 


For  small  crNL,  =  1  —  F, 


so 


PT  =  2  7T  /  (1  —  F)  7(0,  r)  r  dr 

Jo 

POO  POO 

=  2  it  7(0,  r)  r  dr —2  n  F  7(0,  r)  r  dr 
Jo  _  Jo 

‘  V  ^ 

Po 


Substituting  F  back  in, 


(91) 

(92) 


Pt  =  Po  ~  2  tt 


a^NL  7(0,  r)i^^— ^le-/l(0) 


7(0,  r)  r  dr  (93) 


Some  algebra, 


F 


vf  =  1  -  ^  r^NL  72(0,r)Mfl_^Me-M°)  r  dr 

-M3  -Ml  ./n  i-/ 


(94) 


0  Jo 


G 


Reversing  the  approximation,  1  —  G  ~  -yb_, 
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l  -M3  Jo  L 


We  now  turn  to  solving  for  I.  This  is  the  average  intensity.  The  average  intensity 
consists  of  an  J2  term,  which  we  already  defined  above  in  Equation  84.  With  this 
substitution  and  expanding  out  the  definition  for  J2,  we  have 


7=  /Vt4)  [  t  4^^ I  r  dr  (96) 

\P0)  \n2u4(0)J  J0  [J0  u2{L)  L  \  u20 


Combining  the  exponents,  switching  the  order  of  integration  and  careful  algebra 
leads  to  the  following  relation  for  the  average  intensity, 


huj2 (0)  J0  2  +  a2(L-H)2  L 

If  we  let  y 2  =  a2(L  —  77) 2 ,  then  we  can  arrive  at  the  following, 


/  =  _ L  / v  (98) 

Hu2(0)a  L  J^2  +  y2  1  ; 

This  integral  can  now  be  solved  and  we  arrive  at  our  final  form  for  the  average 
intensity, 


Y  Po  V%  ( a  L 

1  =  — t— — — -  arctan  — ■= 
HU2(0)aL[  \2y/2 


Inserting  Equation  99  into  Equation  95, 
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Pt  = 


Pi 


Pi 


1  +  NL  — -  ^ 


VI 

7TUJ2(0 )  a  L 


arcta“  {m 


(100) 


Io 


We  simply  divide  both  sides  by  the  area,  7to;2(0),  to  put  this  in  terms  of  intensity  and 
arrive  at  the  fourth  and  final  intensity  model: 


Method  4  :  I(L ) 


1  +  o’2(f)NL 


Io 


Io 

VI 

a  L 


arctan 


"V" 


(101) 


Vi(a) 

The  average  intensity  accounts  for  the  focused  beam  by  averaging  the  intensity 
profile  to  one  value  throughout  the  length  of  the  test  cell.  Although  this  averaging 
allows  for  an  analytical  form  that  can  be  easily  visualized,  this  model  sorely  under¬ 
estimates  the  absorbance  value  which  in  turn  underestimates  the  TPA  cross-section 
value.  This  is  clearly  shown  in  Figures  33  and  34. 


Comparing  Methods. 

Figures  33  and  34  are  much  like  Figure  16,  but  now  we  include  Methods  3  and 
4.  The  two  figures  share  the  same  trends,  except  for  when  Method  5  crosses  above 
and  then  below  Methods  1  and  2  for  the  a  =  0.555  case  in  Figure  33.  We  see  in 
both  figures  that  Method  3  grossly  overestimates  the  cross-section.  Recall  that  this 
method  only  solves  for  the  intensity  along  the  optical  axis  (/ 3  =  0),  so  it  uses  only 
the  max  value  for  the  intensity  all  along  the  test  cell.  There  is  no  point  at  which 
this  method  is  a  good  approximation.  It  was  surprising  that  this  method  was  so 
divergent  from  the  rest,  although  since  the  second-order  absorption  term  goes  as  I2, 
we  can  understand  where  the  sensitivity  in  I  comes  from.  Method  4  was  the  average 
intensity  model.  Here  we  see  that  the  model  grossly  underestimates  the  TPA  cross- 


77 


section.  For  very  low  absorbance  values  (For  A  less  than  0.03),  this  method  is  a  crude 
yet  acceptable  approximation,  however,  it  quickly  diverges  as  absorbance  increases. 
This  analysis  proves  the  importance  of  having  an  accurate  model  of  the  spatial  profile 
for  the  intensity. 

It  is  also  important  to  note  that  as  a  — >  0,  i.e.  as  the  beam  becomes  more 
collimated,  Methods  1,2,4  &  5  diverge  more  greatly.  Perhaps  the  radial  component 
of  the  intensity  profile  becomes  a  more  significant  factor  in  these  limits. 
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Figure  33.  The  five  methods  compared  against  Method  1,  the  one-photon  Beer’s  Law 
description  of  the  two-photon  absorption  process. 


Table  7  lists  the  eta  factors  for  each  of  the  TPA  models.  Two  focal  lengths  were 
used  in  the  experiment  and  so  two  beam  waists  (to  ~  yielded  two  values  for  a. 

Finally,  Figure  35  is  similar  to  Figure  18  but  includes  the  average  intensity  pre¬ 
diction  for  the  TPA  cross-section.  It  is  important  to  note  that  the  more  complete 
numerical  method  has  less  of  a  statistical  scatter  than  the  average  intensity  method. 
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Figure  34.  The  five  methods  compared  against  Method  1,  the  one-photon  Beer’s  Law 
description  of  the  two-photon  absorption  process. 


This  provided  another  benchmark  pointing  to  the  accuracy  of  our  numerical  method. 
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Table  7.  Comparing  the  77’s  from  the  different  methods. 


Method 

Symbol 

Relation 

a  =  0.555 

a  =  0.303 

2 

772(a) 

772(a) 

1 

1 

3 

'73(a) 

[>LArcTan  (°£)]  (l  +  "2f  ) 

5.88 

3.13 

4 

■/74(a) 

SArcTan  (^) 

0.211 

0.316 

5 

Vb  (r,a) 

*( 

re  “0  (1+“2(z-^)2)  |  ol2L2\  \  1+o2l1) 

J  1 +a2(z-^)2  M  1  4  ' 

See  Figure  14 
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Figure  35.  All  87  cross-section  values  for  both  the  numerical  method  and  the  average 
intensity  method.  It  is  clear  that  the  complete  numerical  method  generates  a  more 
accurate  result.  The  statistical  spread  for  the  average  intensity  model  is  about  twice 
as  large  as  the  numerical  method. 
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Appendix  C.  All  87  Test  Results 


This  study  collected  data  from  87  test  configurations  by  varying  alkali  concentra¬ 
tion,  pump  laser  power  and  beam  waist.  Table  8  displays  all  of  the  doppler-broadened 
cross-sections,  it’s  error  bound,  the  alkali  concentration,  the  pump  power  measured 
just  before  entering  the  optical  beam  train  and  the  beam  waist.  The  right-most  col¬ 
umn  references  the  corresponding  figure  that  displays  the  transformed  data  that  was 
fitted  with  the  Voigt  model  further  below  in  this  appendix. 

Table  8.  All  of  the  87  Doppler-Broadened  Two-Photon  Absorption  Cross-Sections 
Measured  for  the  52Sj/2  — ►  52D5/2  Transition  in  Rb. 


02  (j> 

xlO-21  ^ 

Error 

±  x  10-23  ^ 

N 

N'  1015  a^oms 
cm 3 

Po 

mW 

fim 

Figure 

2.29 

0.56 

1.86 

1200 

69.7 

36 

6.05 

0.66 

3.81 

1200 

69.7 

37 

5.23 

0.52 

3.81 

1230 

69.7 

38 

7.81 

0.99 

4.20 

1230 

69.7 

39 

5.25 

1.02 

7.57 

1230 

69.7 

40 

5.32 

0.99 

7.57 

1230 

69.7 

41 

4.49 

0.45 

1.61 

1200 

69.7 

42 

4.37 

0.64 

1.55 

1200 

69.7 

43 

4.12 

0.73 

1.67 

1200 

69.7 

44 

6.24 

0.89 

2.00 

1200 

69.7 

45 

8.39 

1.40 

2.07 

1200 

69.7 

46 

7.58 

0.91 

2.07 

1200 

69.7 

47 

11.4 

2.89 

2.92 

1200 

69.7 

48 

14.2 

0.84 

3.23 

1200 

69.7 

49 

9.35 

1.14 

3.13 

1200 

69.7 

50 

81 


12.6 

1.07 

4.48 

1200 

69.7 

51 

10.9 

2.58 

4.34 

1200 

69.7 

52 

16.2 

1.31 

4.62 

1200 

69.7 

53 

7.22 

2.09 

5.94 

1200 

69.7 

54 

9.27 

2.26 

5.25 

1200 

69.7 

55 

7.44 

0.65 

5.58 

1200 

69.7 

56 

8.77 

1.02 

5.58 

1200 

69.7 

57 

6.11 

1.04 

2.64 

1200 

69.7 

58 

8.98 

1.17 

3.34 

1200 

69.7 

59 

9.89 

1.13 

4.07 

1200 

69.7 

60 

9.47 

1.56 

4.07 

1170 

69.7 

61 

9.66 

1.51 

4.07 

1150 

69.7 

62 

10.1 

1.31 

4.07 

1150 

69.7 

63 

9.52 

1.34 

4.07 

1130 

69.7 

64 

7.93 

1.51 

4.07 

1110 

69.7 

65 

9.26 

1.36 

4.07 

1110 

69.7 

66 

11.7 

1.55 

4.07 

1090 

69.7 

67 

10.8 

1.74 

4.07 

1090 

69.7 

68 

10.9 

1.59 

4.07 

1070 

69.7 

69 

9.28 

2.22 

4.07 

1070 

69.7 

70 

11.6 

1.45 

4.07 

1050 

69.7 

71 

12.0 

1.54 

4.07 

1050 

69.7 

72 

9.92 

1.35 

4.07 

1030 

69.7 

73 

11.8 

2.58 

4.07 

1030 

69.7 

74 

12.5 

1.26 

4.07 

980 

69.7 

75 

8.99 

1.59 

4.07 

980 

69.7 

76 

82 


13.9 

1.34 

4.07 

930 

69.7 

77 

12.9 

1.40 

4.07 

930 

69.7 

78 

11.8 

1.27 

4.07 

880 

69.7 

79 

10.9 

1.37 

4.07 

880 

69.7 

80 

11.5 

1.56 

4.07 

830 

69.7 

81 

11.7 

1.35 

4.07 

830 

69.7 

82 

10.9 

1.30 

4.07 

780 

69.7 

83 

13.1 

1.54 

4.07 

780 

69.7 

84 

18.4 

2.34 

4.07 

730 

69.7 

85 

15.9 

5.14 

4.07 

730 

69.7 

86 

13.8 

1.59 

4.07 

680 

69.7 

87 

12.1 

1.46 

4.07 

680 

69.7 

88 

14.8 

2.22 

4.07 

630 

69.7 

89 

8.31 

1.35 

4.07 

630 

69.7 

90 

9.19 

1.60 

4.07 

500 

69.7 

91 

9.19 

1.47 

4.07 

500 

69.7 

92 

1.39 

0.31 

2.46 

1200 

94.4 

93 

4.06 

0.41 

2.46 

1200 

94.4 

94 

2.23 

0.34 

2.46 

1200 

94.4 

95 

3.38 

0.42 

3.45 

1200 

94.4 

96 

4.68 

0.40 

3.45 

1200 

94.4 

97 

4.83 

0.57 

4.20 

1200 

94.4 

98 

4.93 

0.73 

4.20 

1200 

94.4 

99 

6.14 

0.63 

4.20 

1200 

94.4 

100 

4.87 

0.49 

6.12 

1200 

94.4 

101 

3.99 

0.76 

6.12 

1200 

94.4 

102 
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5.45 

0.76 

6.12 

1200 

94.4 

103 

5.05 

0.34 

6.12 

1200 

94.4 

104 

4.86 

0.32 

8.52 

1200 

94.4 

105 

3.72 

0.61 

8.03 

1200 

94.4 

106 

2.09 

0.80 

1.29 

1200 

94.4 

107 

6.21 

2.34 

2.55 

1200 

94.4 

108 

4.35 

0.85 

2.55 

1200 

94.4 

109 

3.88 

0.58 

2.55 

1200 

94.4 

110 

4.28 

0.51 

3.12 

1200 

94.4 

111 

4.08 

0.76 

3.12 

1200 

94.4 

112 

5.33 

0.84 

3.34 

1200 

94.4 

113 

7.35 

0.59 

4.20 

1200 

94.4 

114 

6.72 

0.53 

4.20 

1200 

94.4 

115 

6.39 

1.14 

4.20 

1200 

94.4 

116 

5.73 

0.92 

4.20 

1200 

94.4 

117 

5.16 

0.83 

4.20 

1200 

94.4 

118 

7.22 

1.13 

4.20 

1200 

94.4 

119 

8.47 

1.15 

4.20 

1200 

94.4 

120 

10.4 

1.99 

4.20 

1200 

94.4 

121 

9.44 

1.71 

4.20 

1200 

94.4 

122 

The  following  are  the  87  absorption  features  with  transformed  y-axis  and  Voigt 
model  fits.  Please  note  that  the  order  of  the  peaks  is  dependent  on  whether  the  pump 
laser  scanned  up  or  down  in  frequency.  As  an  example,  Figures  36  and  37  have  peak 
#1  (please  refer  back  to  Figure  11)  in  different  locations.  Since  this  work  was  not 
concerned  with  absolute  frequency  position,  no  care  was  taken  to  standardize  the 
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plots  in  this  regard. 


0  2000  4000  6000  8000  10000 

Relative  Frequency,  v,  MHz 


Figure  36.  <t20  =  2.96  x  10"21^. 
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Relative  Frequency,  v,  MHz 


Figure  37.  er20  =  6.05  x  lO"21^. 


85 


0  2000  4000  6000  8000  10000 

Relative  Frequency,  v,  MHz 


Figure  38.  er20  =  5.23  x  10  21 


Figure  39.  =  7.81  x  lO"21^. 
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Figure  40.  cr20  =  5.25  x  10"21^. 
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Figure  41.  cr20  =  5.32  x  10  2l£^- 
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Figure  42.  a2(j)  =  4.49  x  10"21^. 


Figure  43.  <t20  =  4.37  x  10"21^. 


88 


0  2000  4000  6000  8000  10000 

Relative  Frequency,  v,  MHz 


Figure  44.  a24>  =  4.12  x  10"21^. 
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Figure  45.  cr20  =  6.24  x  lO"21^. 


Figure  46.  a2(jl  =  8.39  x  10"21^. 
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Figure  47.  cr20  =  7.58  x  lO"21^. 
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Relative  Frequency,  v,  MHz 


Figure  48.  a2<t>  =  1.14  x  lO"20^. 


Figure  49.  <j2(j>  =  1.42  x  10"20^. 
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Figure  50.  a20  =  9.35  x  lO"21^. 
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Figure  51.  a20  =  1.26  x  lO"20^. 
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0  2000  4000  6000  8000  10000 

Relative  Frequency,  v,  MHz 

Figure  52.  cr2(j>  =  1.09  x  lO"20^. 
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Relative  Frequency,  v,  MHz 

Figure  53.  a2(j>  =  1.62  x  lO"20^. 


93 


0  2000  4000  6000  8000 

Relative  Frequency,  v,  MHz 

Figure  54.  <j2rp  =  7.22  x  10"21^. 
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Figure  55.  cr20  =  9.27  x  lO"21^. 
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Figure  56.  a2l/>  =  7.44  x  10  w 
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Figure  57.  a24>  =  8.77  x  10"21^. 
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Figure  58.  er20  =  6.11  x  10"21^. 
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Figure  59.  cr20  =  8.98  x  10"21^. 
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Figure  60.  cr20  =  9.89  x  lO"21^. 
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Figure  62.  cr20  =  9.66  x  10“21^. 


Figure  63.  a2(j)  =  1.01  x  lO"20^. 
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TPA  Cross-Section,  o2.,  cm  /W 
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Figure  64.  a2(j)  =  9.52  x  10~21^. 


Figure  65.  a2(j)  =  7.93  x  10"21^. 
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Figure  66.  er20  =  9.26  x  lO"21^. 


Figure  67.  cr20  =  1.17  x  lO"20^. 
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TPA  Cross-Section,  a2.,  cm  /W  w  xIO 


Figure  68.  a2<j>  =  1.08  x  lO"20^. 
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Figure  69.  cr20  =  1.09  x  lO"20^. 
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Figure  70.  a2(j>  =  9.28  x  10"21^. 


Figure  71.  cr20  =  1.16  x  lO"20^. 
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0  2000  4000  6000  8000  10000 

Relative  Frequency,  v,  MHz 

Figure  72.  a2<j>  =  1.20  x  lO"20^. 


Figure  73.  a2(P  =  9.92  x  10"21^. 
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Figure  74.  cr20  =  1.18  x  lO"20^. 
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Figure  75.  a2<t>  =  1-25  x  1CT20^. 


104 


0 

2000  4000  6000 

Relative  Frequency,  v,  MHz 

8000 

10000 

Figure  76.  a2<t>  =  8.99  x  10"21^. 

Figure  77.  a2(j)  =  1.39  x  lO"20^. 
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Relative  Frequency,  v,  MHz 

Figure  78.  a20  =  1.29  x  lO"20^. 


Figure  79.  a20  =  1.18  x  lO"20^. 
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Figure  80.  a2(p  =  1.09  x  10"20^. 
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Figure  81.  a2<t,  =  1.15  x  lO"20^. 
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800  ~1  JA 


Figure  82.  a2<j>  =  1.17  x  lO"20^. 
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Figure  83.  cr20  =  1.09  x  lO"20^. 
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Figure  84.  a20  =  1.31  x  lO"20^. 


Figure  85.  a20  =  1.84  x  lO"20^. 
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0  2000  4000  6000  8000  10000 

Relative  Frequency,  v,  MHz 

Figure  86.  a2(/>  =  1.59  x  lO"20^. 
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Figure  87.  a2<j)  =  1.38  x  lO"20^. 
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TPA  Cross-Section,  ct2((),  cm  /W  x-|o‘ 
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Figure  88.  a2(j)  =  1.21  x  lO"20^. 
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Figure  89.  =  1.48  x  lO"20^. 
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Relative  Frequency,  v,  MHz 

Figure  90.  er20  =  8.31  x  10~21^^. 
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Figure  91.  cr20  =  9.19  x  10"21^. 
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Figure  92.  a20  =  9.19  x  10"21^. 
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Figure  93.  a20  =  1.39  x  lO"21^. 
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4  —  21  cm 


Figure  94.  cr2^  =  4.06  x  10  w  . 
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Relative  Frequency,  v,  MHz 


Figure  96.  er20  =  3.38  x  10"21^. 


Relative  Frequency,  v,  MHz 

Figure  97.  er20  =  4.68  x  10~21^. 
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Relative  Frequency,  v,  MHz 

Figure  98.  a2<j>  =  4.83  x  10"21^. 


Figure  99.  a2(j)  =  4.93  x  10"21^. 
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Figure  101.  a2<t>  =  4.87  x  10"21^. 
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Figure  102.  er20  =  3.99  x  10"21^. 


Figure  103.  =  5.45  x  10"21^. 
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Figure  104.  =  5.05  x  10"21^. 


Figure  105.  er20  =  4.86  x  10  21 
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Figure  106.  =  3.72  x  lO"21^. 
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Figure  107.  er20  =  2.09  x  10"21^. 
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Figure  108.  a2<t>  =  6.21  x  lO"21^. 
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Figure  109.  =  4.35  x  10  21 
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Figure  110.  £r2<£  =  3.88  x  10  2l£j%h. 
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Figure  111.  a2</>  =  4.28  x  10"21^. 
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Figure  112.  er20  =  4.08  x  10"21^. 
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Figure  114.  a2(j>  =  7.35  x  lO"21^. 


Figure  115.  a2<t>  =  6.72  x  10"21^. 
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Figure  116.  a2<t>  =  6.39  x  lO"21^. 


Figure  117.  er20  =  5.73  x  lO"21^. 
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Figure  118.  <720  =  5.16  x  10  w 
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Figure  119.  =  7.22  x  10"21^. 
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Figure  120.  a2<t>  =  8.47  x  lO"21^. 
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Figure  121.  a2(t>  =  1.04  x  10"20^. 
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Figure  122.  =  9.44  x  10"21^. 
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Appendix  D.  Voigt  Model  in  IGOR  Pro 


The  Equations  below  are  the  main  pieces  to  the  IGOR  Pro  code  used  for  the  Voigt 
model.  There  are  17  Voigts  corresponding  to  the  17  allowed  hyperhne  transitions. 


’A) 


A  uh  me 2 

cw[2\  = - \j  — - —  ( Lorentzian  width) 


8  kT 


(102) 


cw[  3]  =  — 


mr 


uo  V  2  kT 


(1/ Doppler  width ) 


(103) 


f  Voigt  =  - 


1‘ 


dt 


(104) 


^  J-oo  y  +  (x- 1)2 

In  addition  to  the  equations  above,  the  Voigt  model  contains  the  isotope  fraction, 
the  Boltzmann  factor  and  the  two-photon  relative  signal  strength.  Finally,  the  Voigt 
model  also  includes  four  factors  for  the  3rd  order  polynomial  baseline. 
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